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DEVELOPMENT  AND  EVALUATION  OF  A  MONTE  CARLO  CODE  SYSTEM 
FOR  ANALYSIS  OF  IONIZATION  CHAMBER  RESPONSES 


The  purpose  of  this  work  is  the  development  and  testing  of  a 
Monte  Carlo  code  system  for  calculating  the  response  of  an 
ionization  chamber  to  a  mixed  neutron  and  photon  radiation 
environment.  The  resulting  code  system  entitled  MICAP  -  a  Monte 
Carlo  Ionization  Chamber  Analysis  Package  -  determines  the  neutron, 
photon,  and  total  responses  of  the  ionization  chamber  to  the  mixed 
field  radiation  environment.  The  Monte  Carlo  method  performs 
accurate  simulations  of  the  physical  processes  involved  in  detecting 
radiation  using  ionization  chambers,  and  eliminates  limitations 
inherent  in  approximate  methods. 

The  calculational  scheme  used  in  MICAP  follows  individual 
radiation  particles  incident  on  the  ionization  chamber  wall 
material.  The  incident  neutrons  produce  photons  and  heavy  charged 
particles  and  both  primary  and  secondary  photons  produce  electrons 
and  positrons.  As  these  charged  particles  enter  or  are  produced  in 
the  chamber  cavity  material,  they  lose  energy  and  produce  ion  pairs 
until  their  energy  is  completely  dissipated  or  until  they  escape  the 
cavity.  Ion  recombination  effects  are  included  along  the  path  of 
each  charged  particle  rather  than  applied  as  an  integral  correction 
to  the  final  result.  ENDF/B-V  partial  cross  section  data  have  been 
incorporated  in  the  neutron  transport  module  to  account  for  all 
processes  which  may  contribute  to  the  output  signal.  The  transport 
modules  utilize  continuous  angular  distribution  and  secondary  energy 
distribution  data  when  selecting  the  emergent  direction  and  energy 

x 


of  a  particle.  Furthermore,  reactions  are  treated  as  discrete  and 


allowed  to  occur  with  any  of  the  constituent  nuclides  comprising  a 
mixture.  Finally,  MICAP  incorporates  a  combinatorial  geometry 
package  and  input  cross  section  processors  to  eliminate  restrictions 
in  the  modeling  capability  of  the  code  system  with  respect  to 
geometry,  physical  processes,  nuclear  data,  and  sources. 


To  evaluate  MICAP,  comparisons  were  made  with  results  obtained 
using  other  code  systems  and  with  experimental  results.  Separate 
comparisons  with  other  code  systems  verified  the  validity  of  the 
neutron,  photon,  and  charged  particle  transport  processes  and  the 
nuclear  models  used  to  describe  the  individual  neutron  reactions, 
respectively.  Comparisons  with  mono-energetic  photon  calibration 


a 


experiments  and  with  mixed  neutron  and  photon  radiation  experiments 
verified  the  applicability  of  MICAP  for  analyzing  the  response  of 


ionization  chambers  to  mixed  field  radiation  environments. 
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CHAPTER  I 


INTRODUCTION 

In  the  field  of  radiation  dosimetry,  it  is  often  necessary  to 
establish  an  accurate  relationship  between  a  radiation  field  and  an 
observed  response  in  order  to  infer  physical  quantities  such  as 
radiation  exposure,  energy  transfer,  or  absorbed  dose.  The  rela¬ 
tionship  depends  on  the  characteristics  of  the  radiation,  the  irra¬ 
diated  material,  and  the  detection  device.  Ideally,  the  detector 
should  not  perturb  the  radiation  field,  and  therefore,  provide  an 
observable  response  in  a  known  and  reproducible  manner.  In  prac¬ 
tice,  however,  perturbations  must  be  considered  in  the  interpreta¬ 
tion  of  the  detector  response. 

A  type  of  detector  commonly  used  in  dosimetry  is  the  gas -filled 
ionization  chamber.  The  detector  is  comprised  of  a  container,  a 
gaseous  fill  material,  and  a  charge  collection  system.  The  con¬ 
tainer  and  gas  materials  are  normally  selected  to  closely  match  the 
material  to  be  irradiated  so  as  to  minimize  the  perturbing  effects 
of  the  detector.  Because  of  the  relatively  low  density  gas  region, 
these  detectors  are  referred  to  as  "cavity  ionization  chambers,"  and 
the  methodology  used  to  determine  the  response  of  the  detectors  is 
referred  to  as  "cavity  chamber  theory."^  Using  cavity  chamber 
theory  to  analytically  predict  the  observed  response  of  the  detec¬ 
tors  requires  some  approximate  representations  of  the  physical 
processes  occurring  within  the  detectors. 

The  purpose  of  the  present  work  is  to  develop  and  evaluate  a 
Monte  Carlo  code  system  for  determining  the  response  of  a  gas-filled 


ionization  chamber  in  a  mixed  neutron  and  photon  radiation  environ¬ 
ment.  In  particular,  the  code  system  will  calculate  the  neutron, 
photon,  and  total  responses  of  the  ionization  chamber.  The  Monte 
Carlo  analysis  of  an  ionization  chamber  performs  accurate  simula¬ 
tions  of  the  physical  processes  involved  in  detecting  radiation  and 
eliminates  the  limitations  inherent  in  existing  deterministic 
methods  based  on  cavity  chamber  theory. 

1 . 0  BACKGROUND 

The  cavity  ionization  chamber  is  a  gas -filled  enclosure  in 
which  the  incident  radiation  produces  ionization.  Within  the  enclo¬ 
sure,  there  are  two  or  more  electrodes  which  operate  under  the 
influence  of  an  externally  applied  voltage.  As  the  applied  voltage 
increases,  the  drift  velocities  of  the  electrons  freed  in  the  ioni¬ 
zation  processes  increase  and  ion  recombination  decreases.  At 
saturation  voltage,  ion  recombination  is  at  a  minimum,  yet  the  vol¬ 
tage  is  not  so  strong  as  to  cause  significant  secondary  ionization 
(charge  amplification  through  cascading).  Therefore  the  observed 
signal  is  proportional  to  the  total  energy  deposited  by  charged  par¬ 
ticles  produced  via  the  incident  radiation. 

Most  analytical  and  experimental  techniques  currently  used  in 
radiation  dosimetry  are  based  on  cavity  chamber  theory.  The  funda¬ 
mental  assumption  of  cavity  chamber  theory  is  that  the  dimensions  of 

the  cavity  are  small  compared  with  the  ranges  of  the  electrons  pro- 

2 

duced  in  the  ionization  processes.  More  precisely,  the  theory 
assumes  that  the  size  of  the  cavity  is  such  that: 


1.  the  electron  spectrum  established  in  the  enclosing 


material,  i.e.  in  the  chamber  wall,  is  not  modified  by  the 
presence  of  the  gas  in  the  cavity, 

2.  electrons  generated  in  the  cavity  from  interactions 
between  the  gas  and  the  primary  or  secondary  radiation  are 
negligible,  and 

3.  the  primary  radiation  fluence  (neutron  and/or  photon  flu- 
ence)  is  spatially  uniform  in  the  region  from  which  secon¬ 
dary  electrons  enter  the  cavity. 

The  analytical  and  experimental  techniques  for  photons  (X-rays  and 

gamma- rays)  have  been  extensively  developed  throughout  the  history 
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of  radiation  dosimetry.  Hence,  accurate  determinations  of  the 

absorbed  dose  for  numerous  photon  energies  and  various  source -object 
configurations  are  routinely  accomplished. 

The  determination  of  the  absorbed  dose  associated  with  a  neu¬ 
tron  field  has  not  been  as  extensively  studied  as  that  from  photons. 
High  and  intermediate  "mono-energetic"  neutron  source  experiments 
constitute  most  of  the  work.  Unlike  mono -energetic  photon  sources 
which  emit  photons  at  discrete  energies,  neutron  sources  classified 
as  mono -energetic  actually  emit  a  spectral  distribution  peaked  at 
some  energy.  Therefore,  the  interpretation  of  the  absorbed  dose 
data  requires  knowledge  of  the  incident  neutron  spectrum  and  conse¬ 
quently,  the  accuracy  in  the  determination  of  the  absorbed  dose  may 
be  compromised.  Although  tKe  uses  of  neutron  sources  in  biology  and 
medicine  have  increased  significantly,  comparisons  of  the  results 
from  the  analytical  and  experimental  techniques  used  in  neutron 
dosimetry  continue  to  show  large  discrepancies  in  the  reported 
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absorbed  doses.  These  discrepancies  were  evidenced  at  the  Interna¬ 
tional  Neutron  Dosimetry  Intercomparison  (INDI)^  and  the  European 

g 

Neutron  Dosimetry  Intercomparison  Project  (ENDIP) ,  where  most  of 
the  available  ionization  chambers  (commercial  and  research)  were 
evaluated.  Two  factors  which  contributed  to  the  observed  discrepan¬ 
cies  at  INDI  and  ENDIP  were  the  inconsistencies  in  the  experimental 
procedures  used  to  obtain  the  measured  responses,  and  the  systematic 
differences  in  the  absolute  values  of  the  theoretical  parameters 
used  to  derive  the  doses  from  the  ionization  chamber  measurements . 

The  measurement  of  neutron  dose  is  often  complicated  by  the 
presence  of  a  photon  background  and  because  most  dosimeters  are  sen¬ 


sitive  to  both  neutrons  and  photons. 


Since  the  biological 


effects  of  neutrons  and  photons  are  different,  the  two  components 
must  be  determined  in  order  to  obtain  the  correct  response  of  the 
biological  system  to  the  mixed  field  radiation.  Ionization  chambers 
used  in  dosimetry  work  are  usually  designed  and  constructed  to 
satisfy  the  assumptions  associated  with  cavity  chamber  theory  and  to 
minimize  the  effects  of  approximations  resulting  from  the  theory. 

For  neutrons,  the  maximum  size  of  a  cavity  ionization  chamber  that 
satisfies  the  condition  of  negligible  secondary  particle  production 
is  inconveniently  small.  Consequently,  analyzing  the  data  from 
these  ionization  chambers  using  cavity  chamber  theory  may  produce 
errors  in  the  neutron  dosimetry  results.  Therefore  the  various 
organizations  involved  in  neutron  dosimetry  adopt  different 
approaches  toward  determining  the  absorbed  dose.  These  approaches, 
although  all  based  on  variations  of  cavity  chamber  theory,  have  lead 
to  contradictory  results  and  conclusions. 


‘V. 


1.1  NEED  FOR  THE  PRESENT  WORK 


Recently,  the  Armed  Forces  Radiobiology  Research  Institute 
(AFRRI)  expressed  the  need  for  Monte  Carlo  estimations  of  detector 
responses  in  terms  of  the  electrical  charge  collected  when  the  ioni¬ 
zation  chambers  are  placed  in  mixed  neutron  and  photon  radiation 
environments .  Such  radiation  environments  include  nuclear  battle¬ 
field  environments,  standard  reference  radiation  fields  such  as  that 
of  the  Health  Physics  Research  Reactor  (HPRR)  at  ORNL,  those  at  the 
Army  Pulse  Radiation  Division  (APRD)  and  AFRRI,  and  the  fields  com¬ 
piled  in  IAEA  Report  180,  "Compendium  of  Neutron  Spectra  in  Criti¬ 
cality  Accident  Dosimetry . The  use  of  the  Monte  Carlo  method  is 
generally  regarded  as  the  best  way  to  avoid  the  shortcomings  of  the 
currently  employed  methods  based  on  cavity  chamber  theory. 

Existing  general  Monte  Carlo  code  systems  are  not  tailored  to 
perform  the  ionization  chamber  calculations  for  the  radiation  fields 
of  interest  to  AFRRI.  The  term  "general"  as  applied  to  Monte  Carlo 
indicates  there  are  few,  if  any,  restrictions  in  the  modeling  capa¬ 
bility  of  the  code  system  with  respect  to  geometry,  physical 
processes,  nuclear  data,  and  radiation  sources.  However,  even  a 
general  Monte  Carlo  code  usually  involves  problem- dependent  user- 
written  subroutines  which  tailor  the  code  for  specific  applications. 

There  are  several  instances  in  the  literature  where  specific 

^2 _ *1  ^ 

Monte  Carlo  codes  have  been  written  for  detector  calculations. 

These  codes  are  limited  to  specific  applications  such  as  augmenting 
existing  results  from  analytical  methods,  or  providing  data  needed 
by  the  analytical  methods,  A  code  written  for  a  specific  applica- 


tion  usually  requires  a  major  reprogramming  effort  before  it  is 
applicable  to  another  problem.  Furthermore,  Monte  Carlo  codes  writ¬ 
ten  for  specific  applications  usually  have  the  nuclear  data  (cross 
sections  and  stopping  powers)  included  in  the  program  itself.  This 
further  complicates  the  effort  associated  with  using  the  program  for 
a  different  problem.  The  advantage  of  developing  a  general  Monte 
Carlo  code  system  is  the  adaptability  to  a  wide  range  of  problems 
through  the  use  of  problem- dependent  subroutines.  Also,  general 
Monte  Carlo  code  systems  would  employ  data  pre -processors  to  obtain 
the  nuclear  data  needed  for  a  particular  application. 

1 . 2  PROJECT  OBJECTIVES 

In  summary,  the  objectives  of  the  present  work  were  divided 
into  the  following  tasks: 

1.  Develop  input  data  processors  for  a  neutral  particle  Monte 
Carlo  code  to  arrange  the  pointwise  data  from  the 
Evaluated  Nuclear  Data  File  (ENDF)  library  into  a  format 
compatible  with  the  code. 

2.  Modify  an  existing  neutral  particle  Monte  Carlo  code  to 
calculate  the  physical  processes  occurring  in  a  typical 
ionization  chamber  used  in  mixed  field  dosimetry. 

3.  Develop  a  Monte  Carlo  code  to  calculate  charged  particle 
and  recoil  heavy  ion  energy  loss  processes  in  the  ioniza¬ 
tion  chamber. 

A.  Modify  an  existing  photon-electron  Monte  Carlo  code  to 
calculate  the  physical  processes  occurring  in  the  ioniza¬ 


tion  chamber. 
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5 .  Use  the  new  code  system  to  generate  results  that  can  be 

compared  with  experimental  data  and  with  results  from 

analytical  methods  based  on  cavity  chamber  theory. 

The  PXMORSE^-^  Monte  Carlo  code  (a  continuous  energy  version  of 

MORSE) ^  was  chosen  as  the  code  to  modify  for  the  neutral  particle 
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transport  calculation,  and  the  EGS  photon- electron  transport 

module  of  the  high-energy  calorimeter  Monte  Carlo  code  system  was 
chosen  as  the  code  to  modify  for  the  photon  transport  calculation. 

1 . 3  CALCULATIONAL  PROCEDURE 

The  calculational  scenario  follows  the  individual  radiation 
particles  incident  on  the  ionization  chamber  wall  material.  The 
incident  neutrons  produce  photons  and  charged  particles  (protons, 
alpha  particles,  recoil  ions,  etc.),  and  both  the  primary  and  secon¬ 
dary  photons  produce  electrons.  As  these  charged  particles  and 
electrons  enter  or  are  produced  in  the  chamber  cavity  material  (usu¬ 
ally  a  gas,  but  possibly  a  tissue-like  substance),  they  produce  ion 
pairs  until  their  energy  is  completely  dissipated  or  until  they 
escape  the  cavity.  The  number  of  ion  pairs  produced  is  computed 

using  work  functions,  which  are  defined  as  the  average  energy 
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required  to  produce  an  ion  pair.  The  work  functions  account  for 
charged  particle  and  electron  energy  loss  mechanisms  in  the  detector 
cavity  in  addition  to  ion  production.  The  Monte  Carlo  model 
accounts  for  energy  losses  from  the  cavity  in  the  form  of  delta  rays 
(secondary  low-energy  electrons),  bremsstrahlung  photons,  or  other 
charged  particles  producing  no  ionization.  An  analytical  model  by 


Birks  is  incorporated  to  model  any  ion  recombination  effects  which 
might  occur.  The  Birks  simulation  is  performed  along  each  charged 
particle  and/or  electron  path  rather  than  applied  as  an  integral 
correction  to  the  final  result. 

This  calculational  procedure  models  the  physical  processes 
occurring  in  gas  ionization  more  accurately  than  the  current 
analysis  methods  based  on  cavity  chamber  theory.  All  nuclear  data 
such  as  cross  sections,  stopping  powers,  etc.,  are  utilized  in  a 
pointwise  manner,  which  enhances  the  accuracy  of  the  procedure. 
Continuous  angular  distribution  and  secondary  energy  distribution 
data  are  incorporated  for  selecting  the  emergent  direction  and 
energy  of  a  particle.  Furthermore,  reactions  are  treated  as 
discrete  and  allowed  to  occur  with  any  of  the  constituent  nuclides 
comprising  a  mixture.  Finally,  a  combinatorial  geometry  package  and 
input  cross  section  processors  are  utilized  to  eliminate  restric¬ 
tions  in  the  modeling  capability  of  the  code  system  with  respect  to 
geometry,  physical  processes,  nuclear  data,  and  sources. 

1.4  ORIGINALITY  OF  PRESENT  WORK 

The  final  product  of  the  present  research  is  a  new  code  system 
which  provides  a  unique  capability  in  the  area  of  mixed  field 
dosimetry.  The  present  work  incorporates  the  models  of  all  physical 
processes  occurring  in  gas  ionization  into  the  Monte  Carlo  random 
walk  procedure.  In  particular,  new  and  significant  capabilities  in 
ionization  chamber  response  calculations  are  realized  by  incorporat¬ 
ing  models  describing: 


1.  ion  pair  production  and  charge  collection  processes, 

2.  charged  particle  energy  loss  mechanisms, 

3.  charged  particle  transport  and  ion  recombination  effects, 
and 

4.  nonelastic-nucleus  collisions,  i.e.  (n,p),  (n,d),  (n,t), 
(n,a),  etc. 

Utilizing  all  nuclear  data  such  as  cross  sections,  stopping  powers, 
etc.  in  a  pointwise  manner  yields  a  more  rigorous  treatment  of  the 
particle  transport  processes  than  that  available  from  current 
methods.  By  developing  a  "general"  Monte  Carlo  code  system  tailored 
to  ionization  chamber  calculations,  the  present  work  will  be  useful 
for  many  applications  in  the  field  of  radiation  dosimetry. 

There  is  presently  no  computational  code  system,  Monte  Carlo  or 
otherwise,  that  is  capable  of  performing  the  analysis  of  radiation 
dosimetry  experiments  to  the  extent  developed  in  the  present  work. 
The  "specific"  Monte  Carlo  codes  usually  apply  to  only  one  experi¬ 
mental  set-up,  limit  the  number  of  interactions  allowed  to  occur  in 
the  analysis,  and  restrict  the  experiment  to  simple  geometries,  e.g. 

slab,  concentric  spheres,  or  concentric  cylinders.  The  existing 
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general  purpose  Monte  Carlo  code  systems,  i.e.  MORSE,  MCNP,  and 
23 

TRIPOLI,  have  all  been  written  for  neutral  particle  transport 
analysis.  As  such,  these  codes  concentrate  on  the  interactions 
affecting  the  neutron  and/or  photon  flux  without  regard  to  the  other 
products  (protons,  deuterons,  alpha  particles,  etc.)  of  the  interac¬ 
tions.  Consequently,  these  codes  do  not  model  the  production  and 
transport  of  the  low  energy  charged  particles  essential  to  ioniza¬ 
tion  chamber  response  analysis.  The  present  work  provides  a  new 


computational  tool  which  eliminates  most  of  the  shortcomings  of  the 
methods  currently  being  used  in  ionization  chamber  response 
analysis . 
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CHAPTER  II 


APPLICATION  OF  MONTE  CARLO  TO  THE  SOLUTION 
OF  THE  TRANSPORT  EQUATION 

2.0  BOLTZMANN  TRANSPORT  EQUATION 

The  generalized  time -dependent  integro-differential  form  of  the 
Boltzmann  transport  equation  can  be  derived  by  conserving  particles 
within  a  differential  volume  in  phase  space.  More  specifically,  the 
derivation  equates  the  net  storage  of  particles  within  a  differen¬ 
tial  element  of  phase  space  (drdEdfl)  to  the  particle  gains  minus  the 
particle  losses.  The  derivation  has  been  presented  in  many  texts 
and  publications  and  will  briefly  be  discussed  here.  The  following 

discussion  is  based  on  work  presented  in  S .  N.  Cramer's  disserta- 
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tion  and  in  D.  E.  Bartine's  dissertation. 

The  general  time -dependent  integro-differential  form  of  the 

Boltzmann  transport  equation  is: 

~  -£r  ^(r.E.O.t)  +  n-V<^(r  .E.n.t)  +  Z(r ,  E)rf(r ,  E,0,  t) 

V  ot  t 

( 2  - 1 ) 

-  Q(r,E,0,t)  +  JJdE'dn'Z  (r,E'-*E,n'-»n)^(r,E' ,0' ,t)  , 
where  (r,E,fi,t)  denotes  the  general  seven-dimensional  phase  space, 
r  -  position  variable, 

E  -  the  particle’s  kinetic  energy. 


v  -  the  particle’s  speed  corresponding  to  its  kinetic  energy  E, 


fi  -  a  unit  vector  which  describes  the  particle's  direction  of 
motion, 


t  -  time  variable, 

^(f,E,fi,t)  -  the  time -dependent  angular  flux, 

Q  _  _ 

-  —  <^(r , E,fl, t)dEdfi  -  net  storage  (gains  minus  losses)  per  unit 
volume  and  time  at  the  space  point  r  and  time  t  of  particles 
with  energies  in  dE  about  E  and  with  directions  which  lie  in  dfl 
about  H, 


fl*V^(r , E,0, t)dEdfl  -  net  convection  loss  per  unit  volume  and 
time  at  the  space  point  r  and  time  t  of  particles  with  energies 
in  dE  about  E  and  directions  which  lie  in  dfl  about  fi, 

£t(r,E)  -  the  total  cross  section  at  the  space  point  r  for  par¬ 
ticles  of  energy  E, 

Et(r,E)^(r,E,fl,t)dEdfl  -  collision  loss  per  unit  volume  and  time 
at  the  space  point  r  and  time  t  of  particles  with  energies  in 
dE  about  E  and  directions  which  lie  in  dfl  about  0, 


Zs(r ,E'-*E,Q'-*Q)dEdft  -  the  differential  scattering  cross  section 
which  describes  the  probability  per  unit  path  that  a  particle 
with  an  initial  energy  E'  and  an  initial  direction  fl'  undergoes 
a  scattering  collision  at  r  which  places  it  into  a  direction 
that  lies  in  dfi  about  0  with  a  new  energy  in  dE  about  E, 


2s(r,E'-E,n'-n)<Mr,E'  ,0'  ,t)dE'dfl 


■] 


dEdfl 


inscattering  gain 
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per  unit  volume  and  time  at  the  space  point  r  and  time  t  of 
particles  with  energies  in  dE  about  E  and  directions  which  lie 
in  dfl  about  fl,  and, 


Q(r ,E,fl, t)dEdfl  -  source  particles  emitted  per  unit  volume  and 
time  at  the  space  point  r  and  time  t  with  energies  in  dE  about 
E  and  directions  which  lie  in  dfl  about  0. 

While  Monte  Carlo  codes  are  capable  of  solving  the  time  dependent 
Boltzmann  transport  equation,  the  present  work  is  directed  specifi¬ 
cally  at  the  solution  of  the  static  (time- independent)  Boltzmann 
equation : 


fl.V*(r,E,n)  +  £t(r,E)*(r,E,fl)  -  Q(r,E,fi)  + 


JJdE'dn'Zs(r,E'-E,n'-n)*(r,E' ,«' )  . 


(2-2) 


Equation  2-2  represents  the  most  general  form  of  the  static 
integro-dif ferential  Boltzmann  transport  equation.  As  such,  this 
equation  is  directly  applicable  to  neutron  transport  and  for  most 
applications  involving  photon  transport.  However,  in  ionization 
chamber  response  analyses,  Eq .  2-2  must  be  modified  for  photon  tran¬ 
sport  to  account  for  electron  production  of  bremsstrahlung  photons 
and  photons  produced  via  neutron  interactions.  More  specifically, 
the  inscatter  term  in  Eq .  2-2  becomes: 


ff dE'dfi'2s(r,E'-*E,0'-n)«Hr,E'  ,fi' )  - 
JJdE'dft'Z^(r,E'->E,n'->fi)<^(r,E'  .0')  + 
J/dE'dn'Ze_(^(r,E'-*E,n'-»n)^e(r,E'  ,0'  )  + 


(2-3) 


'uni 


//dE'dO'Sn^(r ,E'-»E,n'-Hl)^n(r ,E'  ,0')  , 

where 

dE'dn,E^(r(E'-EIQ'-n)^^(r,E'  ,0') 
per  unit  volume  at  the  space  point  r  of  photons  with  energies 
in  dE  about  E  and  directions  which  lie  in  dfl  about  0, 

dE' dfl’ Ee^(r  , E' -*E,fl'-*fl)^e(r  ,  E' ,fl' )  j  dEdll  -  bremsstrahlung 
scattering  gain  per  unit  volume  at  the  space  point  r  of  photons 
with  energies  in  dE  about  E  and  directions  which  lie  in  dfl 
about  0, 


dEdfi  -  inscattering  gain 


j^//dE'dn'En^(r,E'-E,n'-n)^n(r,E' ,0' )jdEdn  -  photon  production 
gain  per  unit  volume  at  the  space  point  r  of  photons  with  ener¬ 
gies  in  dE  about  E  and  directions  which  lie  in  dfl  about  fi,  and 


^(r.E.O)  -  the  photon  angular  flux, 


^e(r,E,fl)  -  the  electron  angular  flux,  and, 


<j>  (r.E.O)  -  the  neutron  angular  flux, 
n 

Incorporating  Eq.  2-3  into  Eq .  2-2  and  substituting  <f> ^  for  4>  and  Q 
for  Q  yields: 

O-Vf^r.E.O)  +  Zt(r,E)*  <r,E,0)  -  Q  (r.E.O)  + 

J/dE'dn'E^(r,E'-E,n'-n)^^(r,E’  ,Q’ )  + 

JJdE'dfl'Z  (7,E'-E,0'-0U  (r.E'  ,0’)  + 

G~*y  e 
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J/dE'dO'Z  (r,E'-E,0'-^l)^n(r,E' ,0')  . 
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Equation  2-4  represents  the  generalized  static  integro- differential 
Boltzmann  transport  equation  for  photons  with  bremsstrahlung  produc¬ 
tion. 


As  in  the  case  with  photons,  Eq.  2-2  must  also  be  modified  for 
electron  transport  to  account  for  photon  production  of  electrons  via 
the  photoelectric  absorption  and  pair  production  interactions. 
Furthermore,  at  low  electron  energies,  it  becomes  impractical  to 
simulate  discrete  electron  interactions,  and  Eq.  2-2  must  be  modi¬ 
fied  to  include  a  continuous  energy  loss  term.  This  difficulty  with 
electron  transport  arises  because  the  cross  sections  for  most  elec¬ 
tron  interactions  become  very  large  as  the  electron  energy 
approaches  zero.  The  exact  values  are  not  well  known  and  it  is 
therefore  not  feasible  to  try  to  simulate  every  interaction.  To 
properly  account  for  the  above  considerations,  the  inscatter  term  in 
Eq.  2-2  becomes: 


J/dE'dO'Z  (r,E'-*E,n'-»0)^(r,E'  ,0')  = 

//dE'dn'Se_>e(r,E'-E,n'->0)^e(r,E'  ,fi' )  + 

J/dE'dii'Z  (r,E'->E,n'-n)tf  (r,E',fi')  + 
•y->e  y 


_5 

5E 


[ 


S(EU  (r.E.fl) 
e 


(2-5) 


where 


^JJdE' dfl' Ze_>g(r ,  E'->E  ,0' (r ,  E' ,n' )  j  dEdn  =  inscattering  gain 
per  unit  volume  at  the  space  point  r  of  electrons  with  energies 
in  dE  about  E  and  directions  which  lie  in  dfl  about 


|//dE'dfi'2^e(r,E’-*E,fi'->n)<^(r,E'  ,0' )  j  dEdft  =  photoelectric 
absorption  and  pair  production  gains  per  unit  volume  at  -the 


space  point  r  of  electrons  with  energies  in  dE  about  E  and 
directions  which  lie  in  dfl  about  fl, 

(E)^e(r,E,0)j jdEdfl  -  a  continuous  slowing  down  per  unit 
volume  at  the  space  point  r  of  electrons  with  energies  in  dE 
about  E  and  directions  which  lie  in  dfl  about  0, 

S(E)  -  the  energy  loss  per  unit  pathlength  i.e.,  the  stopping 
power , 

<^e(r,E,fl)  -  the  electron  angular  flux,  and, 

^(r,E,fl)  -  the  photon  angular  flux. 

Incorporating  Eq.  2-5  into  Eq.  2-2  and  substituting  for  4>  and 
for  Q  yields: 

n*^e(r,E,0)  +  2t(r,E)*e(r,E,Q)  -  Qe(r,E,n)  + 

rfdE'dfl'2  (r ,E'-*E,Q'-*Q)^  (r.E’.fi')  + 
e-»e  e 

J/dE’dfl'2  (r.E'-E.fl'-Q)*  (r.E'.O')  + 

7”*e  7 

_  _  1  (2-6) 
S(E)«e(r,E,0)l  . 

Equation  2-6  represents  the  generalized  static  integro-differential 
Boltzmann  transport  equation  for  electrons  with  continuous  slowing 
down  and  photoelectric  absorption  and  pair  production. 

While  Eq .  2-6  is  written  specifically  for  electron  transport, 
it  could  apply  to  any  charged  particle  transport.  The  difficulty  -  r. 
applying  Eq.  2-6  to  other  charged  particles,  i.e.,  protons,  alpha 
particles,  recoil  heavy  ions,  etc.,  is  the  lack  of  cross  section 


data.  The  development  of  theoretical  and  empirical  formulae  to 


simulate  charged  particle  transport  is  not  as  extensive  as  that  for 


electron  transport.  There  are,  however,  numerous  theoretical  formu¬ 
lae  developed  for  calculating  charged  particle  stopping  powers  which 
are  used  in  the  present  work  to  simulate  charged  particle  energy 
deposition  processes.  Using  stopping  powers  to  simulate  the  charged 
particle  energy  deposition  processes  is  acceptable  because  below  20 
MeV  the  probability  for  charged  particle  nuclear  interaction  is 
small  due  to  the  short  tracklengths  of  the  charged  particles.  These 
short  tracklengths  result  from  the  limited  range  of  the  charged  par¬ 
ticles.  Consequently,  the  majority  of  charged  particles  would  depo¬ 
sit  their  energy  before  undergoing  a  collision.  Therefore,  for 
charged  particle  transport,  Eq.  2-6  with  <j>^  substituted  for  <f>  and 

Q  substituted  for  Q  ,  can  be  reduced  to: 
c  e 


Q.V<*c(r,E,fl)  =  Qc(r,E,fl)  +  ^|[s(E)^(r  ,E,0)] 

where 

^c(r,E,n)  -  the  charged  particle  angular  flux,  and, 


(2-7) 


the  definitions  of  the  terms  in  Eq.  2-6  now  apply  for  charged 
particles . 

In  the  Monte  Carlo  method,  a  transport  cross  section  is  used  to 
determine  the  next  collision  site.  Therefore,  a  small  fictitious 
transport  cross  section  is  incorporated  to  force  the  transport  of 
the  charged  particle  to  a  material  boundary.  The  stopping  powers 
are  then  used  to  determine  the  amount  of  energy  deposited  in  that 
material.  This  yields  a  charged  particle  transport  equation  given 
by: 
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(2-8) 


Q.V4  (r.E.O)  +  2*(r ,E)d  (r ,E,fl)  -  Q  (r,E,Q)  + 
c  t  c  c 


//dE'dn'S*(r,E')5(E'^E,fi'^n)^c(r,E' ,0')  + 


aE[s(E)*c(r,E,n)]  , 


where 


*  - 


£t(r,E)^(r,E(fl)dEdfl  -  fictitious  collision  loss  per  unit  volume 
at  the  space  point  r  of  charged  particles  with  energies  in  dE 
about  E  and  directions  which  lie  in  dft  about  fi,  and, 


j^J'JdE’dfl'S*(r,E' )5(E'-»E,n'-n)^c(r,E' ,n' )jdEdfl  -  fictitious 
inscattering  gain  per  unit  volume  at  the  space  point  r  of 
charged  particles  with  energies  in  dE  about  E  and  directions 
which  lie  in  dfl  abour  fi. 

Equations  2-2,  2-4,  and  2-6  represent  the  generalized  static 
integro-differential  Boltzmann  transport  equation  for  neutrons,  pho¬ 
tons,  and  electrons  respectively.  As  such,  these  equations  are 
cross-coupled  because  of  photon  production  of  electrons,  and  neutron 
and  electron  production  of  photons.  Applying  the  Monte  Carlo  method 
to  these  equations  effectively  decouples  the  equations  by  treating 
these  production  terms  as  part  of  the  source  term.  In  other  words, 
Eq .  2-4  can  be  reduced  to  Eq .  2-2  if  the  source  term  is  defined  as: 


QCr.E.O)  -  Q*(r,E,Q)  +  ffdE'dfl'Z  (r , E'-E.O'-fl)*  (r.E'.O')  + 
7  e-*7  e 


J/dE'dO'E  (r.E'-E.O'-^D^Cr.E’  ,0* )  , 


(2-9) 


where  Q^(r,E,ft)  represents  the  external  source  of  photons.  Like¬ 


wise,  Eq .  2-6  can  be  reduced  to  Eq .  2-2  (with  the  added  continuous 
energy  loss  term)  if  the  source  term  is  defined  as: 
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(2-10) 


Q(r,E,0)  -  Q*(r ,E,n)  + 


ffdE'dfl'E  (r,E'-E,0'-n)<£  (r.E',0')  , 
^  7— e  7 


■jSf  —  — 

where  again  Qe(r,E,0)  represents  the  external  source  of  electrons. 


In  light  of  the  above  discussion  and  because  Eq.  2-8  is  a  simplifi¬ 


cation  of  Eq.  2-6,  the  formal  basis  for  the  Monte  Carlo  transport  of 


all  particles  is  provided  by  an  integral  form  of  Eq.  2-2  with  the 


added  continuous  energy  loss  term. 


2.1  RANDOM  WALK  PROCEDURE 


The  previous  section  discussed  the  transport  of  all  particles, 


i.e.,  neutrons,  photons,  electrons,  etc.,  which  contribute  to  the 


ionization  chamber  response  in  a  mixed  field  radiation  environment 


utilizing  the  integro-differential  Boltzmann  transport  equation. 


The  Monte  Carlo  method,  however,  uses  an  integral  form  of  the 


Boltzmann  transport  equation  as  the  formal  basis  for  the  random  walk 


procedure.  A  random  walk  by  a  radiation  particle  is  comprised  of 


its  birth  event,  followed  by  movements  from  one  collision  site  to 


the  next,  and  finally  terminated  by  either  absorption  or  leakage 


from  the  system.  A  reasonable  basis  for  the  Monte  Carlo  random  walk 


procedure  is  the  integral  emergent  particle  density  equation  given 


k  lu 
by  : 


x(r,E,fi)  =  Q(r , E,Q)  + 


Z  (r,E'-*E,0'-*0)  ® 

/JdE'dn'  — - = -  J 


(2-11) 


J dRZ  ( r ,  E ' ) e  ^(r  ,R’ E' ^(r 1  ,  E' , Q' )  , 


2t(r,E')  0 


E 


V 


where 


0<r,R,E,n) 


R 

f  S  (r-R'n.E)dR’ 
0 


is  referred  to  as  the  "optical 


thickness"  and  represents  the  number  of  mean  free  paths  between 
spatial  points  r  and  r' , 


x(r,E,fl)  is  defined  as  the  density  of  particles  leaving  a 
source  or  emerging  from  a  real  collision  with  phase  space  coor¬ 
dinates  (r.E.fl), 

[R]  is  a  spatial  variable  which  relates  a  fixed  point  in  space 
(r)  to  an  arbitrary  point  (r')t  and 

the  definitions  of  the  other  terms  are  given  in  the  discussion 
of  Eq.  2-2  in  the  previous  section. 

The  fundamental  relationship  between  the  emergent  particle  density 
in  Eq.  2-11  and  the  flux  density  is  given  by: 

__  __  __  ____  (2-12) 
x(r,E,0)  -  Q(r,E,0)  +  /JdE'dn'Es(r,E'^E,n'-n)^(r ,E' ,0' )  , 

and  the  flux  density  in  terms  of  the  emergent  particle  density  is 
given  by: 


<j>(r ,  E,ft) 


JdRe^(r’R'E’n) 

0 


*(r,E,n) 


(2-13) 


A  full  discussion  of  the  transformation  of  Eq.  2.2  to  Eq .  2-11  is 
presented  in  Ref.  24  and  will  not  be  repeated  here. 


presented  in  a  simpler  notation  by  defining  the  transport  integral 


operator  as: 


(2-14) 


T(r ' ,R,E,0)  *  JdRZ  (r,E)e 
0  c 


-0(r,R,E,n) 


and  the  collision  integral  operator  as: 


(2-15) 


2  (r,E'-E,n’-n) 
C(r,E'->E,fl'-fi)  -  J/dE'dfi'  — - - - 


St(r,E') 


The  collision  integral  operator  can  be  rewritten  as: 


(2-16) 


2  (r.E'-E.Cl'-ft)  2  (r,E' ) 

C(r,E'-E,n'-n)  -  JJdE' dfl'  — - - -  - 


S.(r,E') 


2t(r,E') 


where 


2s(r,E')  -  J/dEdn2s(r,E'-E,fi'^Q) 


(2-17) 


In  Eq.  2-16,  [2s(r,E'-*E,Q'->fi)/Zg(r  ,E' )  ]  is  a  normalized  joint  proba¬ 


bility  density  function  used  in  selecting  an  emergent  particle's  new 


direction  and  energy  and  [2  (r , E' )/S  (r ,E' ) ]  is  the  nonabsorption 


probability.  Introducing  the  transport  and  collision  operators  into 


Eq .  2-11,  the  integral  emergent  particle  density  equation  in  opera¬ 


tor  notation  is  obtained: 


x(r.E.n)  =  Q(r,E,n)  + 


The  principal  reason  for  selecting  the  integral  emergent  particle 
density  equation  for  the  Monte  Carlo  random  walk  is  because  the 
source  particles  are  introduced  according  to  the  natural  distribu¬ 
tion. 

The  implementation  of  the  random  walk  procedure  is  accomplished 
by  representing  the  emergent  particle  density  x(r,E,Q)  as  a  Neuman 
series : 


x(r,E,0)  -  2  xn(r,E,fl)  , 

n-0 


(2-19) 


where 

*n(r ,E,fl)dEdfl  -  the  density  of  particles  emerging  from  the  nth 
collision  at  the  space  point  r  with  energies  in  dE  about  E  and 
with  directions  which  lie  in  dfl  about  0, 

X^(r,E,fi)  -  the  natural  source  distribution  Q(r,E,n),  and, 

xn(r,E.n)  -  c(^,E'-E,n'-n)T(r'-r,E' ,n')xn_1(^' ,E’ ,0'). 

The  Neuman  series  solution  of  Eq .  2-18  implies  the  following 
sequence  of  events: 

1.  The  random  walk  begins  with  the  selection  of  the 

particle's  phase  space  coordinates,  involving  position 
(Tq),  energy  (Eq) ,  and  direction  (^q) .  according  to  the 
joint  probability  density  function  associated  with  the 


natural  source  distribution  Q(r,E ,0). 


A  flight  distance  R  is  picked  from  the  probability  distri- 

-/SCr.R.Ep.OQ) 

bution  function  (pdF)  Et(r,EQ)e  to  deter¬ 

mine  the  first  collision  site  r^ . 

At  the  collision  site  r^ ,  a  nuclide  from  N  kinds  of 
nuclides  in  the  mixture  is  selected.  The  selection  of  a 
scattering  angle  and  energy  for  a  particular  nuclide  will 
preserve  the  unique  physics  of  each  interaction. 

Once  the  nuclide  has  been  selected,  the  choice  is  made 
between  an  absorption  or  a  scattering  reaction  according 
to  the  nonabsorption  probability  (r^  ,  E^l/E^r^  ,  E^) . 

If  an  absorption  occurs,  the  sequence  is  initiated  again 
for  a  new  particle.  If  a  scattering  reaction  occurs,  a 
new  direction  (0^)  is  selected  according  to  the  marginal 
probability  distribution  function 

JdEEs(r^  ,  Eq-+E  ,0Q->0)/Es(r^  ,  Eq)  .  For  elastic  scattering  and 
inelastic  scattering  to  a  discrete  level,  a  new  energy 
(E^)  is  determined  from  the  kinematic  equations.  For  all 
continuum  reactions,  a  new  energy  (E^)  is  determined  from 
the  conditional  probability  distribution  function 

Repeat  step  2  through  5  until  the  particle  history  is  ter¬ 
minated  due  to  absorption,  escape,  or  because  the 
particle's  emergent  phase  space  coordinates  drop  below 


some  arbitrarily  specified  limit,  i.e.,  energy  cut-off, 
age  cut-off,  etc. 

The  above  sequence  of  events  is  directly  applicable  to  neutron  tran¬ 
sport,  photon  transport,  and  electron  transport  for  electrons  with 
energies  greater  than  an  arbitrary  cutoff  energy  used  to  invoke  the 
continuous  energy  loss  term  shown  in  Eq.  2-6.  However,  for  elec¬ 
trons  with  energies  below  the  cutoff  energy,  and  for  all  charged 
particles,  the  sequence  of  events  changes  to  incorporate  the  con¬ 
tinuous  energy  loss  term.  More  specifically,  the  collision  integral 
operator  defined  in  Eq.  2-15  is  given  by: 

E  _  _  of  )  (2-20) 

_  _  MAX  _  2  (r,E'-E,Q'-fl)  -f  S(E) 

C(r,E,-»E,Q'-*n)  -  f  dE'  /dA'  — - - -  +  2^ - L  , 

I'  St(r,E')  S_(r,E) 

where 

I'  is  an  arbitrary  value  taken  to  be  the  minimum  energy  for 
which  discrete  interaction  is  allowed,  and 

EMAX  is  t^e  max*mum  energy  considered. 

Interactions  involving  particles  with  incident  energy  less  than  I' 
are  considered  to  undergo  a  continuous  slowing  down  with  a  fixed 
energy  loss  per  unit  path  length  traveled.  This  quantity  is 
referred  to  as  the  stopping  power  and  represents  an  integral  over 
the  cross  section  for  the  energy  range  zero  to  I ' .  Therefore,  when 
a  particle  undergoes  a  continuous  slowing  down  type  interaction,  the 
particle  suffers  only  energy  degredation  and  no  angular  deflection, 
i.e.,  straightahead  scattering.  A  complete  derivation  of  the 


J 

I 

I 

continuous  slowing  down  equation  can  be  found  in  Ref.  25  and  will 
not  be  repeated  here. 

Implementing  the  collision  integral  operator  defined  in  Eq.  2-  j 

i 

20  into  the  random  walk  sequence  is  straightforward.  The  sequence  j 

of  events  described  earlier  remains  the  same  for  steps  one  through  I 

four.  Step  five  was  changed  to  include  the  continuous  energy  loss  J 

f 

term  and  the  energy  ranges  over  which  the  two  terms  in  the  collision  j 

integral  operator  apply.  In  particular,  step  five  becomes: 

i 

t 

5.  If  the  incident  energy  of  the  particle  is  greater  than  the  j 

i 

cutoff  energy  I'  and  an  absorption  occurs,  the  sequence  is  ] 

I 

initiated  again  for  a  new  particle.  If  a  scattering  reac¬ 
tion  occurs  and  the  incident  energy  of  the  particle  is 
greater  than  the  cutoff  energy  I'  a  new  direction  (fi^)  is 
selected  according  to  the  marginal  distribution  function 
JdEEg  (r1 ,  EQ-E ,n0-»n)/Ss  ,  Eq  ).  A  new  energy  (E^)  is  deter¬ 

mined  from  the  conditional  probability  distribution  func¬ 
tion  £  (r,  ,  E_-»E ,  )/£  (r-,E_).  If  a  scattering  reaction 

slUUlsiu 

occurs  and  the  incident  energy  of  the  particle  is  less  than 
the  cutoff  energy  I',  the  energy  degredation,  Eq-*E^,  is 
determined  from  the  stopping  power  S(E)  for  the  distance 
transported,  rQ-*r^,  and  the  emergent  direction  is  equal 
to  the  incident  direction  Oq .  If  the  energy  degredation  is 
greater  than  the  incident  energy  of  the  particle,  the  par¬ 
ticle  is  assumed  to  have  slowed  to  rest,  i.e.,  "ranged 


out,"  in  the  medium. 


Step  six  of  the  sequence  remains  the  same.  For  electrons,  the 
cutoff  energy  I'  is  arbitrarily  set  by  the  user  when  processing  the 
cross  sections.  (Typical  values  for  I'  in  the  present  work  were  set 


for  electron  kinetic  energies  of  0.01  MeV)  For  charged  particles 

other  than  electrons,  I'  was  equal  to  Eu.„.  Therefore,  for  these 

MAX 

particles,  all  collisions  were  treated  with  the  continuous  energy 
loss  term. 

An  effect  of  interest  such  as  biological  dose,  energy  deposi¬ 
tion,  or  particle  flux  for  a  given  problem  can  be  expressed  as  any 
one  of  the  following  functionals: 


,  _____  (2-21) 
A  -  /JJP<?(r ,  E ,fi)<£ (r  ,  E ,fl)drdEdfl  , 

,  _____  _  (2-22) 

*  -  //;PV(r,E,n)V-(r,E,n)drdEdn  , 


or 


A  -  JXfPX(r,E,n)x(r,E,n)drdEdn  , 


(2-23) 


where 

6  -  - 

P  (r,E,fl)  -  the  response  function  of  the  effect  of  interest 
due  to  a  unit  angular  flux  (r,E,0), 

P^(r,E,0)  -  the  response  function  of  the  effect  of  interest 
due  to  a  particle  which  experiences  an  event  at  (r,E,Q),  and 

V  —  — 

P  (r,E,fl)  —  the  response  function  of  the  effect  of  interest 
due  to  a  particle  which  emerges  from  a  collision  at  r  with 
phase  space  coordinates  (r,E,n). 


The  cumulative  score  from  a  random  walk  is  the  single  particle 
estimate  of  the  effect  of  interest,  A.  The  Monte  Carlo  estimate  X 
of  the  effect  of  interest  A  is  obtained  by  observing  the  behavior  of 
large  numbers  of  individual  radiation  particles.  More  specifically, 
the  Monte  Carlo  estimate  of  A  is  given  by: 


(2-24) 


5  -  s  • 

i=l 


where  N  is  the  number  of  random  walks  analyzed. 

From  the  Law  of  Large  Numbers,  if  a  true  value  for  A  exists,  then 
the  estimate  X  will  almost  always  approach  A  as  N-*®.  Since  it  is 
impossible  to  transport  an  infinite  number  of  particles,  the  vari¬ 
ance  of  the  Monte  Carlo  estimate  must  be  computed  to  indicate  the 
statistical  uncertainty  of  the  Monte  Carlo  estimate,  X.  The  Monte 
Carlo  estimate  of  the  variance  is  given  by: 


(2-25) 


iPT  ,s ■ 

1=1 


where  N  is  the  number  of  random  walks  analyzed.  It  should  be  noted 
that  as  N-+®,X-*A,  and 


2  LIM  1  1  —  2 

a  =  -  Z  (A.-A) 

N-*®  N  .  ,  l 
1=1 


(2-26) 


In  the  Monte  Carlo  method  as  applied  to  radiation  transport,  the 
standard  deviation  of  the  mean  is  usually  a  reliable  indicator  of 
the  reproducibility  of  X.  The  standard  deviation  of  the  mean  is 
given  by: 


where  N  is  the  number  of  random  walks  analyzed.  From  the  Central 
Limit  Theorem,  there  would  be  a  68%  chance  that  the  X  estimate  of  \ 


rtiques  include  source  biasing,  survival  biasing,  Russian  roulette, 
splitting,  and  exponential  transform.  These  biasing  schemes  have 
been  routinely  employed  in  most  forward  Monte  Carlo  analysis  for 
many  years.  The  primary  function  of  biasing  is  to  encourage  indivi¬ 
dual  random  walks  to  achieve  phase  space  coordinates  from  which  low 
variance  producing  estimates  of  the  effect  of  interest  can  be  made. 
In  the  present  work,  the  basic  programming  for  these  variance  reduc¬ 
tion  techniques  was  included  but  untested.  Therefore,  the  particle 
transport  is  accomplished  using  the  analog  Monte  Carlo  method. 


CHAPTER  III 


NUCLEAR  PROCESSES 

3 . 0  INTRODUCTION 

One  of  the  objectives  stated  in  the  introduction  is  to  develop 
a  general  code  system  that  will  eliminate  most  if  not  all  of  the 
approximations  associated  with  current  methods.  One  avenue  of 
meeting  this  objective  is  to  utilize  all  available  nuclear  data 
applicable  to  the  ionization  chamber  response  analysis. 

3 . 1  NEUTRON  INTERACTIONS 

For  the  neutron  transport  program  this  involves  utilizing  all 
partial  cross  sections,  angular  distributions,  and  secondary  energy 
distributions  currently  available  in  the  Evaluated  Nuclear  Data  File 
-  ENDF/B-V  as  required  for  the  Monte  Carlo  random  walk  procedure. 

The  data  formats  and  procedures  for  processing  the  ENDF/B-V  nuclear 
data  are  presented  in  Ref.  26.  The  recommended  procedures  are 
followed  explicitly  in  the  neutron  cross  section  processors  to 
assimilate  the  data  into  a  format  suitable  for  Monte  Carlo  analysis. 
In  general,  this  involves  formatting  the  cross  sections  into 
linearly  interpolable  cross  section-energy  pairs,  formulating  the 
angular  distributions  into  normalized  energy  dependent  cumulative 
distribution  functions,  and  formulating  the  secondary  energy  data 
into  either  tabulated  probability  distribution  tables  or  into  a  data 
format  suitable  for  sampling  with  one  of  the  ENDF/B-V  secondary 
energy  distribution  functions  (Watt  spectrum,  evaporation  spectrum, 
etc . ) . 
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All  of  the  nuclear  data  in  ENDF/B-V  are  tabulated  for  the 
laboratory  reference  system  except  the  neutron  elastic  and  discrete 
inelastic  scattering  angular  distributions.  These  data  are 
tabulated  for  the  center-of -mass  reference  system  because  these 
neutron  interactions  are  isotropic  in  the  center-of -mass  system  for 
a  wide  range  of  neutron  energies.  Since  the  Monte  Carlo  random  walk 
occurs  in  the  laboratory  system,  exact  angle-energy  relationships 
derived  from  basic  momentum  principles  and  presented  in  Ref.  27  are 
used  to  calculate  the  emergent  neutron  energy  and  direction.  In 
particular,  for  elastic  scattering,  the  scattering  angle  j.n  the 
laboratory  system  is: 

(3-1) 


and  the  emergent  neutron  energy  in  the  laboratory  is: 

!  (3-2) 

E'  ‘  2  +  <1+“>]  ■ 

where 

Mcm  ”  cosine  of  the  scattering  angle  in  the  center-of -mass  sys¬ 
tem  , 

=  cosine  of  the  scattering  angle  in  the  laboratory  system, 

A  =  the  mass  of  the  target  nucleus  divided  by  the  mass  of  a 
neut  ron , 

E  -  the  incident  neutron  energy  in  the  laboratory  system, 


E'  -  the  emergent  neutron  energy  in  the  laboratory  system,  and 


For  discrete  inelastic  scattering,  the  scattering  angle  in  the 
laboratory  system  is  given  by: 


Mlab  “ 


(3-3) 


and  the  emergent  neutron  energy  in  the  laboratory  system  is: 


E' 


E'  + 
cm 


E+2„cn(A+l)(E 

(A+l)2 


(3-4) 


In  Eqs .  3-3  and  3-4,  the  emergent  neutron  energy  in  the  center- of  - 
mass  system  is  needed.  This  can  be  determined  from: 


(3-3) 


where 


E'cm  -  the  emergent  neutron  energy  in  the  center-of-mass  sys¬ 
tem  ,  and , 


Q  -  the  Q  value  of  the  reaction  specified  in  the  ENDF/B-V  data. 
All  other  terms  are  the  same  as  those  used  in  the  elastic  scattering 
equation.  Incorporating  the  exact  angle -energy  formulas  given  by 
Eqs.  3-1  to  3-5  into  the  neutron  transport  program  correctly  models 
elastic  and  discrete  inelastic  neutron  scattering  in  the  laboratory 
system. 


To  account  for  neutron  scattering  at  thermal  energies ,  the  free 
gas  model  is  used  since  it  yields  a  good  approximation  to  the  ther¬ 


mal  flux  spectrum  and  can  be  sampled  without  tables.  The  free  gas 
scattering  model  assumes  that  the  neutrons  are  transported  in  a 
monatomic  gas  having  an  isotropic  Maxwellian  distribution  of  veloci¬ 
ties.  To  obtain  the  emergent  neutron  energy  and  angle,  the  three 


velocity  components  of  the  target  nucleus  are  sampled  from  the 
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Maxwell -Boltzmann  distribution: 


P(V.) 


f— V 

[27rkTj 


exp  - 


AV2 

l 


2kT 


(3-6) 


where 


A  -  the  mass 
neutron, 


of  the  target  nucleus  divided  by  the  mass  of  a 


k  -  the  Boltzmann  constant, 


T  -  the  equilibrium  temperature  (°K)  of  the  target  nuclei,  and 


-  the  target  nucleus  velocity  component  in  the  ith  direction 


for  i  -  x,  y,  or  z. 

The  emergent  neutron  direction  is  sampled  from  an  isotropic  distri¬ 
bution  in  the  center-of -mass  system.  The  emergent  neutron  energy 
and  direction  in  the  laboratory  system  is  then  determined  through 
conservation  of  energy  and  momentum. 

The  emergent  neutron  parameters  for  all  other  reactions  treated 
as  scattering  type  reactions  are  expressed  in  terms  of  post- 
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collision  energy  and  direction  in  the  laboratory  system.  These 
scattering  type  reactions  include  the  continuum  mode  of  inelastic 
scattering,  (n,2n),  (n,3n),  (n, fission),  and  (n,n'x)  where  x 
represents  a  charged  particle  (proton,  alpha  particle,  etc.).  In 
these  reactions,  the  emergent  neutron  angle  is  selected  from  an 
ENDF/B-V  tabulated  angular  distribution  and  the  emergent  neutron 
energy  is  selected  from  an  ENDF/B-V  tabulated  probability  distribu¬ 
tion  or  from  an  evaporation  spectrum  except  for  the  (n, fission) 
reaction  which  uses  an  energy  dependent  Watt  spectrum.  Because  of 
the  structure  of  ENDF/B-V  data,  single  neutron  emission  models  are 
used  for  the  (n,2n),  (n,3n),  and  (n, fission)  reactions  with  the 
weights  of  the  emergent  neutrons  multiplied  by  two,  three  and  »/(E) , 
respectively.  The  parameter,  «/(E) ,  is  the  average  number  of  neu¬ 
trons  produced  per  fission  event  tabulated  as  a  function  of  energy. 

All  charged  particle  production  via  absorption  reactions  i.e., 
(n,p),  (n,a),  (n,d),  etc.,  or  inelastic  scattering  reactions  with 
charged  particle  emission  assume  isotropic  emission  of  the  charged 
particle  in  the  laboratory  system.  The  charged  particle  emergent 
energy  is  selected  from  a  general  evaporation  spectrum  given  by: 28 


F(E-E' ) 


E/ 
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-E'/e 
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(3-7) 


where  I  is  a  normalization  constant  dependent  on  the  nuclear  tem¬ 
perature  ( 6 )  and  the  emergent  particle  energy  range  (E'  .  ,E'  ) 

min  max 

and  6  is  given  by: 28 


e  -  4.016  x  10- 


max 
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(3-8) 
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with  E  -  E  +  Q  -  CB  in  units  of  eV.  In  Eqs.  3-7  and  3-8,  E  is 
max 

the  incident  neutron  energy,  Q  is  the  ENDF/B-V  Q  value  for  the  reac¬ 
tion,  CB  is  the  Coulomb  barrier,  A  is  the  mass  number  in  atomic  mass 

units,  and  E'  is  the  charged  particle  exit  energy  (with  the  Coulomb 

29 

barrier  added  back  in).  The  Coulomb  barrier  is  calculated  from: 


(3-9) 


CB  -  0.75 


Zl  Z2  e 

r*  r 
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where 


Z ^  -  the  charge  number  of  the  recoil  heavy  ion, 


Z^  -  the  charge  number  of  the  charged  particle, 


e  -  the  charge  of  the  proton,  and 


r^  and  x^  -  the  atomic  radii  for  the  recoil  heavy  ion  and 
charged  particle  respectively. 

The  atomic  radius  r^  as  a  function  of  mass  number  is  given  by  the 
following  table: 


2  <  A  <  4 


1.20  x  10" 


2.02  x  10 


2.43  x  10" 
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1.70  A. 1/3  x  10~13 

l 


•If 

$ 

I 


•!a 

•!k 


I 


$ 

$ 


I 


The  factor  0.75  is  set  arbitrarily  to  account  for  charged  particle 
emission  below  the  Coulomb  barrier. 

All  compound  nuclei  excited  through  neutron  interactions 
(except  elastic  scattering)  possess  the  potential  of  emitting  one  or 
more  photons  while  decaying  to  ground  state.  Therefore,  for  every 
neutron  interaction  except  elastic  scattering,  a  check  is  made  to 
determine  if  secondary  photon  production  data  are  available.  If 
there  are  data  available,  a  photon  is  produced  with  a  direction 
assumed  isotropic  in  the  laboratory  system,  an  energy  selected  from 
the  ENDF/B-V  tabulated  distribution,  and  weight  equal  to  the 
ENDF/B-V  energy  dependent  multiplicity.  There  are  some  materials 
with  absorption  reactions  resulting  in  ground  state  transitions. 

For  these  reactions  a  test  is  implemented  to  insure  no  secondary 
photon  generation  occurs. 

The  photon  production  capability  is  programmed  to  model  the 
natural  physical  processes  as  accurately  as  possible.  Unfor¬ 
tunately,  the  photon  production  data  in  ENDF/B-V  are  not  well  known 
for  some  of  the  individual  neutron  interactions  in  many  materials. 
Consequently,  for  these  materials,  the  ENDF/B-V  photon  production 
data  have  been  lumped  into  one  file  encompassing  the  individual  neu¬ 
tron  interactions  which  might  produce  a  secondary  photon.  The  neu¬ 
tron  transport  code  utilizes  these  data  to  produce  a  photon  whose 
direction,  energy,  and  weight  are  representative  average  values  for 
these  neutron  interactions. 

Each  neutron  interaction  produces  a  recoil  heavy  ion.  In  ioni¬ 
zation  chamber  response  analysis,  recoil  heavy  ions  can  deposit 
energy  in  the  active  medium  thereby  contributing  to  the  detector 


signal.  Therefore,  the  energy  and  direction  of  the  recoil  heavy 
ions  must  be  determined  for  each  neutron  interaction.  This  is 
accomplished  using  energy  and  momentum  balances  for  all  incident  and 
exit  particles  produced  in  the  collision. 

3 . 2  PHOTON  INTERACTIONS 

The  photoelectric  effect,  Compton  scattering  and  pair  produc¬ 
tion  are  the  three  processes  considered  in  the  photon  transport  pro¬ 
gram.  In  the  photoelectric  process,  a  photon  is  absorbed  and  an 
orbital  electron  is  emitted  with  approximately  all  of  the  energy  of 
the  photon  transferred  to  the  electron.  The  electron  is  ejected 
from  its  shell  with  kinetic  energy  equal  to  the  incident  photon 
energy  less  the  binding  energy  of  the  orbital  electron  (with  a  very 
small  amount  of  energy  taken  up  by  the  nucleus) .  Momentum  is  con¬ 
served  by  recoil  of  the  whole  atom.  The  photoelectric  process  is 
the  dominant  photon  interaction  for  photons  with  energies  below  0.1 
MeV.  Furthermore,  because  of  momentum  and  energy  conservation  con¬ 
siderations,  the  photoelectric  process  is  more  likely  to  occur  with 
the  more  tightly  bound  K- shell  and  L- shell  electrons. 

Compton  scattering  is  the  elastic  scattering  of  a  photon  by  an 
essentially  free  electron.  The  incident  photon  imparts  energy  to 
the  electron  and  emerges  from  the  collision  with  a  new  direction  and 
energy.  In  the  Compton  scattering  process,  the  electrons  are 
assumed  to  be  free  -  that  is,  not  bound  within  the  atom  or  interact¬ 
ing  among  themselves.  The  Compton  effect  becomes  important  for 
incident  photon  energies  greater  than  0.1  MeV. 


In  the  pair  production  process,  a  photon  interacts  with  the 
electrostatic  field  of  a  charged  particle.  The  incident  photon  is 
completely  absorbed  resulting  in  the  production  of  an  electron- 
positron  pair.  This  photon  interaction  requires  an  incident  photon 
energy  of  at  least  twice  the  rest  mass  of  an  electron,  i.e.,  1.022 
MeV,  and  becomes  increasingly  more  probable  as  the  photon  energy 
increases  above  1.022  MeV.  Momentum  and  energy  is  conserved  in  the 
process  by  the  recoil  of  the  charged  particle.  Either  the  nucleus 
or  an  atomic  electron  can  provide  the  necessary  electrostatic  field 
for  the  process.  Because  of  the  greater  charge  of  the  nucleus,  the 
interactions  are  primarily  nuclear.  If,  however,  the  process  occurs 
in  the  field  of  an  atomic  electron,  energy  is  transferred  to  the 
recoil  electron  and  a  triplet  is  produced.^  The  triplet  consists 
of  the  recoil  electron  and  the  electron-positron  pair.  Because  of 
momentum  and  energy  considerations,  the  threshold  energy  for  triplet 
production  is  four  times  the  rest  mass  of  an  electron,  i.e.,  2.044 
MeV. 

The  interaction  cross  sections  for  the  photoelectric,  Compton, 
and  pair  production  processes  exhibit  a  Z  dependence .  The  Compton 
effect  is  linearly  dependent  on  the  number  of  electrons  per  atom  and 
therefore  is  proportional  to  Z.  The  photoelectric  process  is  pro¬ 
portional  to  Zn,  where  n  varies  from  3  to  5,  and  diminishes  rapidly 
with  increasing  photon  energy  for  any  element.  The  pair  production 

process  is  proportional  to  Z  for  interaction  with  atomic  electrons 

2 

and  proportional  to  Z  for  interactions  with  the  nucleus.  However, 
only  in  iron  and  higher  Z  materials  does  the  pair  production  process 

account  for  a  significant  portion  of  the  energy  absorption  for 
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incident  photon  energies  below  10  MeV.  Consequently,  most  ioniza¬ 
tion  chamber  response  analysis  for  photons  will  be  dominated  by  the 
photoelectric  and  Compton  processes.  The  theoretical  and  empirical 
formulae  used  to  generate  the  cross  sections  for  the  above  photon 
interactions  are  presented  in  detail  in  Ref.  19. 

3 . 3  ELECTRON  INTERACTIONS 

Elastic  scattering  off  the  nucleus,  inelastic  scattering  off 
the  atomic  electrons,  bremsstrahlung  production,  and  positron 
annihilation  are  the  electron  (and  positron)  interactions  considered 
in  the  photon  transport  program.  As  stated  earlier,  difficulties 
arise  in  electron  transport  because  the  cross  sections  for  all  the 
above  interactions  (except  annihilation)  become  very  large  as  the 
electron  energy  approaches  zero.  Therefore,  it  is  not  feasible  to 
simulate  every  interaction,  and  they  are  lumped  together  and  treated 
in  a  continuous  manner.  Cutoff  energies  are  used  to  distinguish 
between  continuous  and  discrete  interactions.  Any  electron  interac¬ 
tion  which  produces  a  delta- ray  with  an  energy  greater  than  the 
electron  cutoff  energy  or  a  photon  with  an  energy  greater  than  the 
photon  cutoff  energy  is  considered  a  discrete  event.  All  other 
electron  interactions  are  treated  in  a  continuous  manner  giving  rise 
to  continuous  energy  losses  and  direction  changes  to  the  electron 
between  discrete  interactions. 

The  continuous  energy  loss  models  complicate  electron  transport 
because  the  cross  section  varies  along  the  path  of  the  electron  and 
the  electron  path  is  not  straight.  This  complication  is  rectified 
by  introducing  an  additional  fictitious  interaction  which  forces  the 


total  cross  section  to  remain  constant  along  the  electron  path.  If 
this  fictitious  interaction  is  chosen,  the  electron  experiences  a 
straight-ahead  scattering  i.e.,  no  interaction  at  all.  Otherwise, 
the  interaction  is  considered  real  and  is  dealt  with  in  the  usual 
manner . 

Finally,  to  account  for  multiple  scattering  of  the  electrons, 
the  transport  of  the  electron  between  interactions  is  divided  into 
small  steps.  Within  each  small  step  the  electron  is  assumed  to  fol¬ 
low  a  straight  line,  and  multiple  scattering  is  accounted  for  by 
changing  the  electron's  direction  at  the  end  of  each  small  step. 

The  size  of  the  steps  is  kept  small  enough  to  insure  that  the  true 
electron  path  length  is  not  much  larger  than  the  straight  line  path 
length,  and  that  the  angle  between  the  initial  and  final  directions 
follows  the  appropriate  angular  distribution. 

The  above  discussion  is  condensed  from  material  presented  in 
Ref.  19.  As  in  the  case  with  photon  interactions,  the  theoretical 
and  empirical  formulae  used  to  generate  the  cross  sections  for  the 
above  electron  interactions  are  presented  in  detail  in  Ref.  19. 

3.4  CHARGED  PARTICLE  INTERACTIONS 

In  solving  the  Boltzmann  transport  equation  for  charged  parti¬ 
cles,  i.e.,  protons,  alpha  particles,  recoil  heavy  ions,  etc.,  the 
only  interactions  considered  in  the  present  work  are  continuous 
energy  loss  mechanisms.  This  is  due  to  the  lack  of  charged  particle 
cross  section  data  for  most  materials.  To  simulate  a  Monte  Carlo 
random  walk,  a  very  small  fictitious  transport  cross  section  is 
incorporated  to  determine  the  next  "collision"  site.  Because  this 


cross  section  is  very  small  (1.0  x  10'^-®),  the  charged  particle 
transport  kernel  will  almost  always  provide  the  distance  to  the 
material  boundary  as  the  distance  to  be  transported.  Prior  to  tran¬ 
sport,  the  particle's  range  for  complete  energy  loss  is  calculated 
from  the  stopping  power  and  range/energ_  table  using  the  following 
relation: ^ 


R(E)  -  / 


E  dE 
o  dE/dx  ’ 


(3-10) 


where 

R(E)  -  the  range  of  a  charged  particle  in  a  particular  medium 
with  energy  dE  about  E,  and, 

dE 

—  -  the  stopping  power  of  a  particular  medium  for  a  charged 
particle  with  energy  dE  about  E. 

If  the  particle's  range  is  less  than  the  distance  to  the  material 
boundary,  all  the  particle's  energy  is  assumed  deposited  in  the 
material.  If  the  range  is  greater  than  the  distance  to  the  material 
boundary,  the  stopping  power  and  range/energy  table  are  used  to 
determine  the  fraction  of  energy  deposited  in  the  material  and  the 
particle  is  transported  to  the  material  boundary  with  the  remaining 
energy.  This  process  is  repeated  until  the  particle  loses  all  of 
its  energy.  The  theoretical  formulae  used  in  the  present  work  for 
generating  the  stopping  powers  and  range/energy  tables  are  given  in 
Ref.  32. 

Recombination  effects  are  sometimes  a  problem  in  gas  or  liquid 
ionization  chambers.  Recombination  results  from  ionized  electrons 


recombining  with  charged  ions  thereby  decreasing  the  charge 


collected.  If  recombination  effects  are  significant,  the  nonlinear¬ 
ity  of  the  charge  collected  (Q)  is  taken  into  account  using  Birks' 


(3-11) 


dx  1+kB^ 

dx 


where  kB  is  the  recombination  constant  determined  experimentally. 

If  kB  is  zero,  the  charge  collected  is  proportional  to  the  energy 
deposited,  i.e.,  In  the  code  system,  recombination  effects 

are  computed  for  each  individual  particle  when  they  are  determined 
to  be  significant.  However,  in  ionization  chamber  analysis,  recom¬ 
bination  effects  are  usually  negligible.  This  formulation  for  cal¬ 
culating  recombination  effects  was  chosen  because  of  the  ease  of 
application  to  the  stopping  powers  and  range/energy  tables  used  in 
the  charged  particle  transport. 


CHAPTER  IV 


SYSTEM  VERIFICATION  AND  RESULTS 

4.0  INTRODUCTION 

Accomplishment  of  the  first  four  objectives  outlined  in  Chapter 
I  led  to  the  development  of  MICAP  -  A  Monte  Carlo  Ionization  Chamber 
Analysis  Package.  MICAP  is  a  modular  code  system  designed  to  calcu¬ 
late  the  response  of  an  ionization  chamber  to  a  mixed  neutron  and 
photon  radiation  environment.  Figure  4-1  depicts  a  flow  diagram  of 
the  overall  calculational  procedure  employed  in  MICAP.  For  this 
collection  of  codes,  the  approach  has  been  to  develop  a  "general" 
ionization  chamber  analysis  code  system,  applicable  to  a  wide  range 
of  problems,  through  the  use  of  input  arrays  and  user  supplied  sub¬ 
routines.  A  user's  manual  for  MICAP  is  given  in  Ref.  33  complete 
with  sample  problems  for  each  module. 

Establishing  the  validity  of  the  nuclear  models  incorporated 
into  the  transport  modules  of  the  code  system  involved  implementa¬ 
tion  of  a  five  part  verification  program.  The  verification  program 
included  comparing  results  obtained  with  existing  code  systems  and 
comparing  with  experimental  results.  Each  part  of  the  verification 
program  investigated  a  particular  aspect  of  the  radiation  transport 
processes  applicable  to  ionization  chamber  analysis.  More  specifi¬ 
cally,  the  verification  program  involved: 

1.  Comparisons  with  MORSE^  results  to  determine  validity  of 
neutron  transport,  secondary  photon  production,  and  photon 


transport  processes. 


2.  Comparisons  with  MACK-IV^  and  RECOIL^  results  to  deter¬ 
mine  validity  of  nuclear  models  used  to  describe  indivi¬ 
dual  neutron  reactions. 

3.  Comparisons  with  05S-^  results  to  determine  validity  of 
charged  particle  transport  processes. 

4.  Comparisons  with  experimental  results  for  photon  calibra¬ 
tion  experiments. 

5.  Comparisons  with  experimental  results  for  mixed  neutron 
and  photon  radiation  experiments . 

The  remainder  of  this  chapter  discusses  each  part  of  the  verifica¬ 
tion  program. 

4.1  COMPARISON  WITH  MORSE 

To  verify  the  neutron  transport,  secondary  photon  production, 
and  photon  transport  processes  in  MICAP,  an  iron  slab  transmission 
problem  was  analyzed  with  both  MORSE  and  MICAP.  The  MORSE  code  sys¬ 
tem  was  chosen  for  the  comparison  because  its  performance  has 
already  been  verified  through  comparisons  with  experiments  and  com¬ 
parisons  with  other  code  systems . , 17 , 37  Consequently,  comparisons 
with  MORSE  should  provide  a  reliable  indication  of  the  validity  of 
the  neutron  transport,  secondary  photon  production,  and  photon  tran¬ 
sport  processes  in  MICAP  The  calculations  involved  modeling  a  neu¬ 
tron  point  source  incident  on  the  front  face  of  iron  slabs  having 
thicknesses  10  and  20  cm  The  material  parameters  for  the  iron  slab 
are  presented  in  Table  4-1  and  the  geometrv  configuration  is  shown 
in  Figure  4-2  Both  code  systems  calculated  the  neutron  and 


Table  4-1.  Material  Parameters  for  the  Iron  Slab 


Parameter 

Iron 

Chemical  formula 
Density 

Iron  Atom  density 


Value 

Fe 

3 

7.870  gm/cm 

8.490  x  10‘2  at/b-cm 


ORNL-DWG  86-15408 


Figure  4-2.  Geometry  Configuration  for  the  Iron  Slab. 
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secondary  photon  leakage  spectra  out  the  back  face  of  the  iron  slab, 


and  the  secondary  photon  production  in  the  iron  slab. 


The  iron  cross  section  data  used  in  MICAP  differed  from  those 


used  in  the  MORSE  calculations.  MICAP  used  ENDF/B-V  point  cross 


section  data  thinned  to  within  a  tolerance  of  0.1%.  The  MORSE  cal¬ 


culations  for  the  10  cm  iron  slab  used  an  ENDF/B-IV  105n-2l7  cou¬ 


pled,  P^,  multigroup  library.  °  To  insure  different  cross  section 


data  bases  were  not  a  source  of  discrepancy,  the  MORSE  calculations 


for  the  20  cm  iron  slab  used  an  ENDF/B-V  37n-21y  coupled,  P^ ,  multi¬ 


group  library^®  as  well  as  the  ENDF/B-IV  105n-2l7  coupled,  P^ ,  mul- 


tigroup  library. 


The  10-cm-thick  iron  slab  was  analyzed  using  a  mono-directional 


point  source  centered  on  the  front  face  of  the  slab  and  the  20-cm- 


thick  iron  slab  was  analyzed  using  an  isotropic  point  source  cen¬ 


tered  on  the  front  face.  Both  iron  slabs  were  analyzed  using  a  U 


fission  spectrum  neutron  source  and  a  mono -energetic  neutron  source 


(14.2  MeV  for  the  10  cm  slab  and  14.0  MeV  for  the  20  cm  slab).  The 


MORSE  calculations  used  Russian  roulette  and  splitting  in  the  tran¬ 


sport  and  the  MICAP  calculations  were  performed  using  analog  Monte 


Carlo.  All  calculations  used  100,000  source  neutrons  incident  on 


the  iron  slab. 


The  neutron  and  secondary  photon  leakage  flux  spectra  out  the 


back  face  of  the  10  cm  iron  slab  and  the  secondary  photon  production 


spectra  are  presented  in  Figures  4-3  to  4-5  respectively  for  the 


U  fission  neutron  source.  The  figures  indicate  excellent  agree¬ 


ment  between  MICAP  and  MORSE  for  the  back  face  leakage  spectra  and 


secondary  photon  production  spectra.  The  spectral  differences  seen 
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Figure  4-5.  Comparison  of  MICAP  and  MORSE  Secondary  Photon 
Production  Spectra  in  a  50  cm  x  50  cm  x  10  cm  Iron  Slab  with  a 
Point  Mono-Directional  Fission  Neutron  Source. 


in  Figures  4-3  to  4-5  are  within  the  statistical  errors  quoted  for 
the  two  calculations.  The  statistical  errors  for  the  two  calcula¬ 
tions  are  unacceptable  for  benchmarking  purposes,  however,  they  do 
provide  a  reasonable  comparison  of  the  transport  processes  in  MICAP 
and  MORSE.  Similar  results  are  shown  in  the  comparison  of  the  MICAP 
and  MORSE  neutron  and  secondary  photon  leakage  flux  spectra  out  the 
back  face  of  the  20  cm  iron  slab  and  the  secondary  photon  production 
spectra  presented  in  Figures  4-6  to  4-8  for  the  isotropic  14.0  MeV 
mono-energetic  neutron  source.  While  the  results  again  show  spec¬ 
tral  differences  between  MICA?  and  MORSE,  comparisons  of  the 
integral  quantities  presented  in  Tables  4-2  and  4-3  for  the  two  10 
cm  iron  slab  cases  and  Tables  4-4  and  4-5  for  the  two  20  cm  iron 
slab  cases  show  excellent  agreement.  In  Table  4-6,  the  integral 
quantities  for  the  MORSE  calculations  using  the  ENDF/B-IV  and 


i 


ENDF/B-V  cross  section  data  sets  are  compared  with  the  MICAP 
ENDF/B-V  results  for  the  20  cm  iron  slab  and  the  isotropic  14.0  MeV 
mono-energetic  neutron  source.  The  fission  neutron  source  results 
were  not  compared  because  of  significant  differences  between  the  two 
input  fission  source  spectra.  The  only  significant  difference 
between  the  ENDF/B-IV  and  ENDF/B-V  MORSE  results  in  Table  4-6,  is 
the  lower  secondary  photon  production  calculated  in  the  MORSE 
ENDF/B-IV  results.  The  lower  secondary  photon  production  conse¬ 
quently  yields  a  lower  integral  photon  leakage  flux  and  current  out 
the  back  face  of  the  slab.  While  the  integral  quantities  in  Table 
4-6  show  slight  differences  between  the  MORSE  ENDF/B-IV  results  and 
the  MICAP  and  MORSE  ENDF/B-V  results,  the  spectral  results  (not 

i 

‘  shown)  were  comparable  and  within  the  statistical  errors  of  the 
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Figure  4-6.  Comparison  of  MICAP  and  MORSE  Neutron  Leakage 
Flux  Spectra  out  the  Back  Face  of  a  50  cm  x  50  cm  x  20  cm  Iron 
Slab  with  a  Point  Isotropic  14.0  MeV  Neutron  Source. 
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Table  4-2.  Comparison  of  MICAP  and  HORSE  Integral  Quantities 
for  the  50  cm  x  50  cm  x  10  cm  Iron  Slab  with  a  Point 
Mono-Directional  Fission  Neutron  Source 


Quantity 

MICAP 

MORSE 

Neutron  flux  leakage  out  the  back  face  (n/cm^) 

3.15-04 

(0.007)3 

3.07-04 

(0.007) 

Neutron  current  leakage  out  the  back  face  (n/cm^) 

1.97-04 

(0.004) 

1.90-04 

(0.004) 

Secondary  photon  production  (7/ n) 

7.85-01 

(0.013) 

7.89-01 

(0.004) 

Photon  flux  leakage  out  the  back  face  (y/cm^) 

5.38-05 

(0.046) 

5.30-05 

(0.020) 

Photon  current  leakage  out  the  back  face  ( 7/cm 2) 

3.03-05 

(0.052) 

2.96-05 

(0.014) 

cl  —  — 

Fractional  standard  deviation  (<r(x)/x). 


Table  4-3.  Comparison  of  MICAP  and  MORSE  Integral  Quantities 
for  the  50  cm  x  50  cm  x  10  cm  Iron  Slab  with  a  Point 
Mono-Directional  14.2  MeV  Neutron  Source 


Quantity 

MICAP 

MORSE 

Neutron  flux  leakage  out  the  back  face  (n/cm^) 

3.92-04 
(0. 007)a 

3.78-04 

(0.007) 

Neutron  current  leakage  out  the  back  face  (n/cm^) 

2.50-04 

(0.004) 

2.47-04 

(0.004) 

Secondary  photon  production  (7/n) 

2.46+00 

(0.004) 

2.34+00 

(0.003) 

Photon  flux  leakage  out  the  back  face  (y/cm^) 

2 . 18-04 
(0.015) 

1.90-04 

(0.014) 

Photon  current  leakage  out  the  back  face  (y/cm^) 

1.21-04 

(0.011) 

1.06-04 

(0.010) 

ei  —  — 

Fractional  standard  deviation  (a(x)/x) . 


Table  4-4.  Comparison  of  MICAP  and  MORSE  Integral  Quantities 
for  the  50  cm  x  50  cm  x  20  cm  Iron  Slab  with  a 
Point  Isotropic  Fission  Neutron  Source 


Quantity 

MICAP 

MORSE 

2 

Neutron  flux  leakage  out  the  back  face  (n/cm  ) 

1 . 08-04 
(0 . 013)a 

9.05-05 

(0.014) 

2 

Neutron  current  leakage  out  the  back  face  (n/cm  ) 

6.43-05 

(0.008) 

5.30-05 

(0.005) 

Secondary  photon  production  (y/n) 

5.57-01 

(0.006) 

6.14-01 

(0.005) 

2 

Photon  flux  leakage  out  the  back  face  (y/cm  ) 

4.33-06 

(0.058) 

4.80-06 

(0.046) 

2 

Photon  current  leakage  out  the  back  face  (y/cm  ) 

2.59-06 

(0.051) 

2.84-06 

(0.040) 

fractional  standard  deviation  (er(x)/x)  . 


Table  4-5.  Comparison  of  MICAP  and  MORSE  Integral  Quantities 
for  the  50  cm  x  50  cm  x  20  cm  Iron  Slab  with  a 
Point  Isotropic  14.0  MeV  Neutron  Source 


Quantity 

MICAP 

MORSE 

2 

Neutron  flux  leakage  out  the  back  face  (n/cm  ) 

1.49-04 
(0 . 017)a 

1.43-04 

(0.012) 

2 

Neutron  current  leakage  out  the  back  face  (n/cm  ) 

8.80-05 

(0.008) 

8.56-05 

(0.007) 

Secondary  photon  production  (y/n) 

3.36+00 

(0.003) 

3.43+00 

(0.003) 

2 

Photon  flux  leakage  out  the  back  face  (y/cm  ) 

3.23-05 

(0.026) 

3.38-05 

(0.027) 

2 

Photon  current  leakage  out  the  back  face  (y/cm  ) 

1.97-05 

(0.022) 

2 .05-05 
(0.022) 

Table  4-6.  Comparison  of  MICAP  and  MORSE  Integral  Quantities 
for  the  50  cm  x  50  cm  x  20  cm  Iron  Slab  with  a  Point 
Isotropic  14.0  MeV  Neutron  Source  and  using  the 
ENDF/B-IV  and  ENDF/B-V  Cross  Section  Data  Bases 


y 

.« 
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calculations.  Therefore,  the  different  cross  section  data  bases 
were  not  considered  significantly  different  for  the  purposes  of  the 
comparison  study.  The  final  verification,  offered  in  Figure  4-9, 
compares  the  PHOTON  module  and  MORSE  secondary  photon  flux  leakage 
spectra  out  the  back  face  of  the  20  cm  iron  slab  for  the  isotropic 
14.0  mono-energeti c  neutron  source  distribution.  In  this  com¬ 
parison,  the  secondary  photon  production  source  generated  in  MORSE 
was  transported  by  both  MORSE  and  the  PHOTON  module  of  MICAP.  Once 
again,  Figure  4-9  shows  that  the  spectral  differences  are  within  the 
statistical  error.  These  differences  were  evidenced  in  all  four 
calculations;  however,  comparisons  of  the  integral  quantities  in 
Table  4-7  show  good  agreement  for  all  cases.  Based  on  these  com¬ 
parisons  and  the  fact  that  MORSE  is  considered  a  viable  code, 
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Figure  4-9.  Comparison  of  PHOTON  and  MORSE  Photon  Leakage 
Flux  Spectra  out  the  Back  Face  of  a  50  cm  x  50  cm  x  20  cm  Iron 
Slab  with  a  Point  Isotropic  14.0  MeV  Neutron  Source  using  the 
MORSE  Generated  Secondary  Photon  Production  Data. 


Table  4-7.  Comparison  of  PHOTON  and  MORSE  Integral  Photon 
Leakage  out  the  Back  Face  of  the  Iron  Slab  Using  the 
MORSE  Generated  Secondary  Photon  Production  Data 
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the  neutron  transport,  secondary  photon  production,  and  photon  tran¬ 
sport  processes  in  MICAP  are  considered  valid 


4.2  COMPARISON  WITH  MACK- IV  AND  RECOIL 

The  neutron  and  photon  leakage  flux  spectra  and  secondarv  pho 
ton  production  spectra  presented  in  the  last  section  demonstrated 
the  validity  of  the  radiation  transport  processes  The  comparisons 
however,  did  not  verify  all  aspects  of  the  nuclear  models  used  to 
describe  the  individual  neutron  reactions  important  to  ion  i  ::a*  i  or. 
chamber  response  analysis.  More  specifically,  verification  is 
needed  for  the  nuclear  models  describing  the  energy  transfer  front 
the  neutron  to  the  reaction  products  for  all  partial  cross  sections 
used  in  MICAP.  To  accomplish  this  goal,  a  problem  consisting  of  . 
mono-directional  14.2  MeV  mono-energetic  neutron  source  ir,<  idet  • 
a  40-cm-thick  iron  slab  was  calculated  using  MI'"AP  The  , 

parameters  and  geometry  configuration  were  the  same  as  •ha’ 
the  spectral  comparisons.  The  MICAP  calculation  us.-d  ui,  . 
neutrons  and  set  an  energy  cut-off  at  1*  Me';  •<«  nt.lv  ,t ;  „ 

collisions.  The  40-cm- thickness  of  iron  >  s 1  •*»-.» r  f  :  .  •  ■ 

insured  an  adequate  number  of  source  int.-r-c  •  •  .  . ■ 

charged  particles  and  recoil  heavy  ions  ;  :  n  -  .  •: 
module  were  processed  in  the  HEAVY  module 
kinetic  energies  of  the  charged  part  i-  l*-s  • 
duced  by  the  individual  neutron  reactin'  •  :  • 

energies  calculated  in  program  HEAVY  we : . 
appropriate  cross  sections  (at  1  .  Me' 


for  the  individual  neutron  reactions.  These  average  kinetic  ener¬ 
gies  and  kerma  factors  were  then  compared  to  results  obtained  with 


MACK- IV34  and  RECOIL.35 

The  kerma  factor  comparisons  involve  using  microscopic  cross 
section  data  evaluate'  at  14.2  MeV.  Because  MACK- IV  and  RECOIL  use 
ENDF/B-IV  cross  section  data  and  MICAP  uses  ENDF/B-V  cross  section 
data,  the  microscopic  cross  sections  for  the  individual  neutron 
reactions  used  in  all  three  programs  are  presented  in  Table  4-8  As 

seen  in  Table  4-8,  there  are  considerable  differences  in  some  of  the 
cross  sections.  The  cross  section  differences  between  MACK- IV  and 
MICAP  are  attributed  to  differences  in  the  cross  section  data  bases 
ENDF/B-IV  versus  ENDF/B-V)  and  to  the  energy  weighting  models  used 
•-o  create  the  MACK- IV  multigroup  cross  section  set  The  differences 
between  MACK  IV  and  RECOIL  are  primarily  attributed  to  differences 
it  cross  section  processing  models 

The  results  of  the  average  kinetic  energv  and  kerma  factor  cal- 
ul.it  ions  performed  in  MICAP  show  good  agreement  with  results  from 
MACK  IV  and  RECOIL  The  comparison  with  MACK  IV  Table  *  rj  i  lumps 
•he  leroi!  nucleus  and  <harged  particle  results  for  a  paificulat 
•*-<'•  i  'ii  ’  nge  •  h  e  r  whe  reas  "he  comparison  with  RECOI  L  Table  «  1 ; '  1 

•  ■  -  ,  jde  s  ’ t  e  'barged  pa  r '  :  ■  1  e  r  e  s  i  I  *  s  entirely  In  viewing  fables 

•  *fd  *  *  he  <-  f  f  e  ’  s  of  ♦  lie  r  o«.s  se-  *  1  of  da*  -1  b.»se  .»  t  •-  t  • 
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Table  4-8.  Iron  Microscopic  Cross  Sections  (at  14.2  MeV)  Used 
in  the  MICAP,  MACK- IV,  and  RECOIL  Kerma  Factor  Calculations 


1  ross  Sec  t  ion  (  b ) 

RECOIL 

MACK- IV 

MICAP 

T  fal 

2.54+00 

2  .  57-t-OO 

1 

. 5  7+00 

V. 1  a  s  t  i  c 

1 .  16+00 

1 . 17+00 

1 

. 2  3+00 

i  i  ■ .  2  n ) 

4.57-01 

4.32-01 

4 

. 7  5-01 

m  . n'a) 

3.67-03 

1.99-03 

2 

.  88-0  5 

1  >  . n'p) 

4.13-02 

3.40-02 

4 

.05-02 

( n , n ' )  resolved 

1.16-01 

1 .17-01 

l 

.11-01 

i  :  , n’ )  continuum 

6.01-01 

6.32-01 

'■> 

.  f'6-01 

■  •  y ) 

2.08-04 

<.  .  10-04 

8 

.  1  0—0  . 

■ 1  •  P) 

1.01-01 

1.25-01 

1 

2  2-01 

-i  .  d ) 

1.90-02 

1.89-02 

6 

.  / 4-0  5 

'  i  • .  r ) 

4.80-04 

1 . 34-04 

6 

58-04 

ii'.  He ) 

3  .  50-04 

1 . 30-04 

6 

38-0 

n  ,  f>  ) 


3.92-02 


3  98-02 


5  .  9/, -02 


Table  4-9.  Comparison  of  MICAP  and  MACK- IV  Average  Kinetic  Energies 
(E)  and  Kerma  Factors  for  the  Partial  Reactions  at  14.2  MeV 
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Table_4-10 .  Comparison  of  MICAP  and  RECOIL  Average  Kinetic  Energies 
(E)  and  Kerma  Factors  for  the  Partial  Reactions  at  14.2  KeV 


RECOIL 

MICAP 

RECOIL 

MICAP 

Kerma  Factor 

Kerma  Factor 

Reaction 

E  (eV) 

E  (eV) 

(b-eV) 

(b-eV) 

Total 

2.24+05 

2.28+05 
(0 . 002)b 

5.69+05 

5.86+05 

Elastic 

1.09+05 

1.08+05 

(0.003) 

1.26+05 

1 . 33+05 

(n, 2n) 

2.68+05 

2.82+05 

(0.006) 

1.23+05 

1 . 19+05 

(n,n'ar)a 

4.28+05 

5.92+05 

(0.069) 

1 . 57+03 

1.70+03 

(n,n' p)a 

2.69+05 

3.19+05 

(0.019) 

1 . 11+04 

1.29+04 

( n . n ' )  resolved 

3.66+05 

3.63+05 

(0.011) 

4 .25+04 

4.04+04 

f  n , n ' )  continuum 

3.31+05 

3.39+05 

(0.005) 

1 . 99+05 

2 . 02  +  04 

in,  i ) 

2.52+05 

2.52+05 

(0.117) 

5.24+01 

2.04+07 

a 

(n,p) 

3.42+05 

3.40+05 

(0.010) 

3.45+04 

4  1  J+0.4 

, ,a 

1  n  ,  d  ) 

3  .  55+05 

4.03+05 

(0.040) 

6 .  75+03 

7.71+ni 

a 

i  n  .  t  ) 

3.28+05 

3 .93  +  05 
(0.25  ) 

1 . 57+02 

7  40  +  0(1 

3U  a 

1  n  ,  He ) 

3  88+05 

4.05+05 

1 . 16+07 

7  1 8  •  <  d 

(0.240) 

a 

i  ti.a) 

6  .  31+05 

8 . 49+04 

7  •  /  +  ()<» 

•'  1  .  + "  . 

(0.018) 

Fx i t  charge 

part ic le 

not  i nc 1 uded 

Fractional  standard  deviation  (o(x)/>:  i 


1 


results  are  seen  In  Table  4-10.  In  general,  the  nuclear  models  used 
in  MICAP  to  describe  the  individual  neutron  reactions  are  adequate 


for  calculating  ionization  chamber  responses.  The  principal  differ¬ 
ences  noted  above  are  due  to  different  cross  section  data  bases  and 
cross  section  processing  techniques. 

4.3  COMPARISON  WITH  05S 

In  the  last  section,  verification  of  the  nuclear  models  used  to 
describe  the  individual  neutron  reactions  was  described  In 
essence,  this  procedure  verified  the  charged  particle  and  recoil 
heavy  ion  source  generation  processes  in  program  NEUTRON  in  much  the 
same  way  as  the  comparisons  to  MORSE  verified  the  secondary  photon 
production  processes  The  charged  particle  and  recoil  heavy  ion 
transport  processes  in  program  HEAVY,  however,  remain  to  he  veri 
fied.  To  accomplish  this  task,  a  comparison  to  OSS ^  was  made  ('is 
is  a  specialized  Monte  Carlo  code  for  calculating  pulse  h<  i gh t  lis 
trihutions  due  to  mono  -  energet  i  e  neutrons  incident  on  organii  m  it! 
tillators  f)‘>S  is  specialized  in  that  the  ■  ode  handles  only 
carbon  hydrogen  organic  compounds  and  models  only  i igh*  ■  it  i!  it 
■  V  1  i  tide  r  s  U‘>s  is  f  ur  f  tie  1  lest  lifted  ill  ’lie  •  a  t  hot.  nut  t:  vd  i  opt  ' 
neutron  react  ions  which  ate  allowed  ’  o  oi.'it  t  <  t  '  r  ;i  « !  - 1 

test!  jet  i  oils  .  O  i  S  has  beet:  used  ex  t  etis  i  ve  1  •,  •  . .  i  .  ■  •  i  !  . '  •  ;  i  . 

he  i  gh t  d  i  s  t  i  i  tni '  i  otis  t  o  i  .  a  1  i  h  i  t  i  ot  .  -  pe  i  i  in  '  1  „  1  •  !  .... 

t  ie  u  t  I  otts  A  I  I  lint  >  a  I  i  son  of  pulse  heigh'  d  i  s  •  i  ;  1  .■ ,  ■  ,  ■  •  . 

1  hi  S  and  M  1 1  A 1'  would  t  he  i  ••  f  o  t  i  i  if  ’hi  fit  »  •  •  *  .  i  • 
heavy  ion  rianspot’  models  in  piopi  np  HEAV'i 


n; 


The  calculations  involved  modeling  a  mono-directional  14.2  MeV 


mono-energetic  neutron  source  distributed  over  the  face  of  a  polvvi 
nyltoluene  (grade  BC501)  plastic  scintillator  The  material  parame 
ters  for  BCbOl  are  presented  in  Table  4-11  and  the  source  and 
geometry  configuration  are  presented  in  Figure  4-10  OSS  is  dif¬ 
ferent  from  typical  Monte  Carlo  programs  in  that  the  number  of 
source  neutrons  undergoing  at  least  one  collision  is  input  rather 
than  the  total  number  of  source  particles  tracked  OSS  then  tracks 
source  particles  until  the  required  number  have  undergone  at  least 
one  collision  To  obtain  adequate  statistics  for  the  pulsi  height 
distribution  comparison  tfsd's  less  than  ♦  10*>  approximately 
40  (t()()  individual  source  particles  unde  t  going  at  least  one  «  nl  1  ;  %  i«.t 
were  considered  sufficient  Therefore  bo  tit  ')'>S  and  MI»AP  processed 
approximately  1  '0.000  incident  soutce  neut . ons 

Tiie  comparison  of  t  lie  pulse  height  distributions  gene  t  at  «  d  it 


O'S  and  MI<AP  is  presented  graphically  i  n  figure 


Both  <ll  s'  ! 


'''H  Inns  ale  normalized  *  o  bit  OOO  <  <u  •  t  ii  bu t  i  ng  snu  1  •  i-  |u!  r  !.  .  (  -.  ,  >  ■<* 
plot  ted  is  turn  !  inns  of  energy  in  oba  1  *  units  •  if*  o t ■  1  •  mi' 

he  base  line  i  n  t  e  r  '  e  )■  •  "  t  I  s  •  t  •  t  gll  *  i  Ut  f  1  ’  '  .  '  1  a  I  ;  ,  .  ■  ,<  * 
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Table  4  11  Material  Parameters  tor  t  lie  Rr  Ml ; 
Pol  yv  iny  1 1  oluene  Plastic  Si  1  nr  i  1  !  n  t  or 


Par  a  me  *  e  • 


.  •  l  r -  .  ,  ’  • *  ■  it 
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POLYVIMYLTOLUENE 


Figure  4-10.  Geometry  Configuration  for  the  BC501 
Polyvinyltoluene  Pleetlc  Scintillator. 
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Figure  4  11  Coaperlaon  of  NICAP  and  05S  Pulse  Heigh' 
Distributions  for  a  BC501  Polyvinyltoiuene  Plastic  Sclntlila*or 
with  a  Mono  Directional  14  ?  NeV  Neutron  Source 


Figure  4-11  shows  excellent  agreement  between  05S  and  MICAP.  Three 
distinct  contributions  can  be  seen  in  the  pulse  height  distributions 
presented  in  Figure  4-11,  with  each  contribution  dominated  by  a  dif¬ 
ferent  particle  type.  The  low  energy  contribution  (below  4  x  10'^ 
cobalts)  is  due  primarily  to  carbon  recoil,  the  intermediate  energy 
contribution  (between  4  x  10" 2  and  approximately  8  x  10" *  cobalts) 
is  due  primarily  to  alpha  particle  recoil,  and  the  high  energy  con¬ 
tribution  (above  approximately  8  x  10"!  cobalts)  is  due  primarily  to 
proton  recoil.  These  contributions  were  plotted  separately  in  Fig¬ 
ures  4-12  to  4-14  to  show  better  resolution.  The  carbon  recoil 
(Figure  4-12)  and  proton  recoil  pulse  heights  (Figure  4-14)  agree  to 
within  the  statistics  of  the  Monte  Carlo  calculations.  The  alpha 
particle  recoil  pulse  heights  (Figure  4-13)  however,  reflect  differ¬ 
ences  in  the  models  employed  in  the  two  programs 

The  two  neutron  reactions  generating  the  alpha  particle  recoil 
pulse  height  are  the  ^C(n.n'3o)  and  the  ^C(n.a)  reactions  In 
viewing  Figure  4-1  1 ,  the  pulse  height  for  OSS  contains  two  distinct 
peaks  at  approximately  1.0  x  10" ^  and  2.S  x  10' '  cobalts  and  a  small 
peak  at  approximately  /  .  S  x  10"  ^  cobalts  The  two  large  peaks  are 
produced  bv  alpha  particles  generated  from  the  ' -Y  n.ti'lc*'  reaction 
and  the  small  peak  is  produced  bv  alpha  particles  generated  t  roir  tin 
'■<  in. u)  reaction  The  two  large  distinct  peaks  result  from  the  O  ->  S 
•leatment  of  ^-Ytn.n'lai  collisions  Carbon  contains  eighteen 
discrete  levels  of  inelastic  scattering  capable  of  decaying  to 
g i  mind  state  via  three  alpha  particle  emission  To  simplify  pro¬ 
gramming  tnd  save  core  storage  space.  OSS  allows  only  the  three  most 
probable  lev*  Is  t  <>  occur  Two  of  these  levels  .  -11  '>'•<*  of  tin-  t  i"< 


LIGHT  ENERGY  (cobolts) 


Figure  4-12.  Comparison  of  MICAP  and  OSS  Pulse  Height 
Distributions  due  Predominantly  to  Carbon  Recoil  for  a  BC501 
Polyvinyl toluene  Plastic  Scintillator  with  a  Mono- Directional 
14.2  MeV  Neutron  Source. 
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Figure  4-13.  Comparison  of  MICAP  and  05S  Pulse  Height 
Distributions  due  Predominantly  to  Alpha  Recoil  for  a  BC501 
Polyvinyltoluene  Plastic  Scintillator  with  a  Mono-Directional 
14.2  MeV  Neutron  Source. 
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Figure  4-14.  Comparison  of  MICAP  and  05S  Pulse  Height 
Distributions  due  Predominantly  to  Proton  Recoil  for  a  BC501 
Polyvinyltoluene  Plastic  Scintillator  with  a  Mono-Directional 
14.2  MeV  Neutron  Source. 

(in  05S)  and  correspond  to  the  two  distinct  peaks  shown  in  Figure 
4-13.  MICAP  maintains  all  discrete  levels  of  inelastic  scattering 
data  and  therefore  generates  much  broader  peaks  in  the  alpha  parti 
cle  pulse  height  distribution  as  seen  in  Figure  4-13.  The  MICAP 
pulse  height  distribution  is  considered  a  more  accurate  represents 
rion  of  the  true  pulse  height  distribution.  While  the  spectral 
shapes  in  the  alpha  particle  pulse  height  are  different,  the 
integral  contributions  are  in  excellent  agreement 

Additional  comparisons  between  03S  and  MICAP  were  made  on  the 
zero  bias  efficiency  (defined  as  the  number  of  source  particles 
undergoing  at  least  one  collision  divided  bv  the  total  number  oi 
source  particles)  and  on  the  average  number  of  collisions  per  col¬ 
liding  neutron.  05S  calculated  a  zero  bias  efficiency  of  0.3.'  and 
1.41  collisions/colliding  neutron.  The  MICAP  results  were  in  exce 


lent  agreement  with  a  zero  bias  efficiency  of  0.36  and  1.46 
collisions/colliding  neutrons.  Even  though  MICAP  required  more  exe¬ 
cution  time,  (6  minutes  cpu  for  MICAP  versus  1  minute  cpu  for  05S  on 
an  IBM  3033)  the  more  detailed  nuclear  models  and  the  non- restricted 
application  capability  of  MICAP  make  it  a  viable  detector  analysis 
code  system. 

4.4  COMPARISON  WITH  PHOTON  CALIBRATION  EXPERIMENTS 


The  preceding  comparisons  with  existing  codes  have  demonstrated 
the  various  components  of  the  radiation  transport  processes  impor¬ 
tant  in  ionization  chamber  response  analysis.  To  demonstrate  the 
ability  of  using  the  general  Monte  Carlo  code  system  MICAP  to  calcu¬ 
late  ionization  chamber  response  functions,  two  mono - energe t i c  pho 
ton  calibration  experiments  were  analyzed.  The  first  experiment 
calculated  was  the  Health  and  Safety  Laboratory  tHASL)  calibration 
experiment  of  the  high  pressure  argon  ionization  chamber  systems 
used  in  the  measurement  of  environmental  radiation  exposure  rates 
The  material  properties  for  the  HASL  high  pressure  argon  chambet  at* 
presented  in  Table  4-12  and  the  geometry  model  is  presented  in  t i , 
ure  4-15.  The  second  experiment  calculated  was  'he  Armed  r  •>:  «  ■ 
Radiobiology  Research  Institute  t  A  ERR  1  >  <>  .  *  1  :  •  ■ :  *  •  •  .  .  • 

ionization  chambers  performed  at  AFRK!  and.  .  H  •  , 

Standards  t  NBS  >  *'"  *•  The  AFRK  I  ioni/at  i>»t 


graph  i  t  e/carbon  dioxide  ■  ■  >  «rd  'issu«  e;  ..  .  .  -  • 

equivalent  gas  i  TK  TK  1  ■>0  cm’  chambe  •  s  the  t...  •  .  -■ 
and  ( TE /TE )  0  5  cm ^  chambers  and  *  h«  Tl 


Table  4-12.  Material  Parameters  for  the  Health  and 
Safety  Laboratory  (HASL)  High  Pressure  Argon 
Ionization  Chamber 


Parameter 

Aluminum 

Chemical  formula 
Density 

Aluj., ;  num  atom  density 

Porce lain 

Chemical  formula 
Dens i  C  v 

Aluminum  atom  densitv 
Ox v gen  atom  der.sitv 
Si  1 i o on  atom  de  ns i t  v 

Stainless  Steel  30-* 
Chemical  formula 
Dens i t  v 


Value 


Al 

2  6  9  Q  g.cm^ 

6.02'  x  10  ■  at  bom 

Al  'Oi  .  Si'1' 

*-  > 

2  t<00  ►tre  :: 


ORNL-DWG  86-15403 


ALL  DIMENSIONS  ARE  0  797  - - H 

IN  CENTIMETERS  .0  890^ 


Figure  x  lh  Gfomftrv  Configuration  for  the  Health  and 
Uh'ratorv  HASI.'  High  Pressure  Argon  Ionization  Chamber. 
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in  Table  4-13  and  the  geometrv  models  *  • 

0.05  cm^  chambers  are  shown  it  K  i 
tively.  Both  sets  of  experiments  p!,i.  •  ■>* 

collimated  beam  of  mono  -  ene  rge  t  !■  ;(;>'*•>•  .• 

charge  collected  (which  is  propm  '  i  ••••-« ! 
the  gas).  The  experiments  were  p.,s  :  •  ;  e<: 

effects  from  the  experimental  rooms  we  >  .  .  , 

The  calculations  modeled  '},>  ;■>.  .  -  ,• 

Figures  4-15  to  4-18  exact  lv  us  i  t:g  ■  .  • 

input  included  the  geomet  t  v  de,<  ; 
subprograms.  The  first  subprogram  to.*.  . .  ■ 
but  ion  generated  in  program  PHOTPRFP  »  •  - 
energy,  and  direction.  The  second  su(q  i  ,• 
transport  program  PHOTON,  oal-ula'ed  » t  • 
sited  in  the  cavitv  hv  each  mntribir  ;  -  ;  , 

'he  results  to  obtain  the  to*  t ;  .-nerg  .  ••  ;  • 

it:  is  quant  i  t  v  was  then  norma  i  ;  r.>  d  :  •  •  •  , 

cd  compared  to  the  experimental  "t.  . 

■ '  •  i  pt:>  i  •  on  .  ross  section  da '  a  we  t  .  g.  •  .  •  , 
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Table  A- 13.  (continued) 


Parameter 

Luc i te 

Chemical  formula 
Density 

Hydrogen  atom  density 
Carbon  atom  density 
Oxygen  atom  density 

Micarta 

Chemical  formula 
Dens i ty 

Hydrogen  atom  density 
Carbon  atom  density 
Oxygen  atom  density 

Carbon  Dioxide  @  1  Atm 
Chemical  formula 
Dens i r  y 

Carbon  atom  density 
•ixvgeri  atom  density 

\  1  r  («i  1  Atm 

Chemical  formula 
Dens i t  v 

Nitrogen  atom  density 
Oxvgen  atom  density 
Argon  atom  density 

[issue  Krjui valent  Gas  (3  1  Atm 
1  hemic a  1  formula 
[tens  i  t  v 

Mvdrogpn  atom  density 
arhon  atom  density 
Ni'rogeii  atom  density 
'  .  gen  a  1  i  ■  m  dens  i  t  y 

■  n-  r ,  j. ; ;  i  .•  n  f  I’  1  a  s  t  i .  A  1)0 

i  i  f  "  t  mu  1  a 


■  *  •  1 !  i  •  orr,  dens  i  '  v 

•  :  '  •  "  i  '  '  >IT,  dele,  |  V 

'■  ’  ’■  O’  <  •  t  '  ■  ,!!■  de  !  S  ;  • 

g  e  •  i  •  ■  nr  1*  1 1  s  1  i  . 


Value 


C5H8O2 

1 . 200  gm/cm^ 

5.777  x  10'^  at/b-cm 
3.611  x  10'2  at/b-cm 
1.444  x  10~2  at/b-cm 


C12H10O2 
1 . 300  gm/cm^ 

4.206  x  10"2  at/b-cm 
5.047  x  10'2  at/b-cm 
8.412  x  10'3  at/b-cm 


C02 

1.965  x  IQ'3 


gm/cmj 


2.690  x  10"3  at/b-cm 
5.380  x  10‘3  at/b-cm 


0.756  N  +  0.232  0  +  0.013  Ar 
1.292  x  10' 3  gm/cm3 
4.203  x  10'3  at/b-cm 
1.128  x  10‘3  at/b-cm 
2.503  x  10' 7  at/b-cm 


0.644  CH4  +  0.324  C02  +  0.032  N2 

1.138  x  10"3  gm/cm3 

6.937  x  10 '3  at/b-cm 

2.602  x  10 " 3  at/b-cm 

1.713  x  10 '6  at/b-cm 

1  .  7o4  x  10' 3  at/b-cm 


0.102H  *  0  ,68  C  +  0.059  0  4 
0  036  N  *  0.018  Ca  4  0.017  F 
!  12  7  gin  cm  ‘ 

6  H  '  1  v  ] o  ■  at  /b -  cm 
t  142  x  1  o  •’  at  b-cm 

1  x  In  '  at  b  cm 

2  ><)..  >.  lo  !  at  b  cm 

6  n‘f,  v  !  o  *  ,,t  h  cm 

i  0,»u  .;  ;  !  I  *  at  b  cm 


Is 


Figure  4-16.  Geometry  Configuration  for  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  50  cm3  Ionization  Chamber. 
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Figure  4-18.  Geometry  Configuration  for  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  0.05  cm  Ionization 
Chamber . 
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Table  4-14  Comparison  of  Measured  and  Calculated  Detector 
Responses  for  the  Photon  Calibration  Experiments  of  the 
Health  and  Safety  Laboratory  (HASL)  High  Pressure 
Argon  Ionisation  Chamber 


Phot  on 

Source 

Source 

Energy 

(MeV) 

Detector 

Pressure 

( Atm ) 

Calculated 
Response 
( A//iR  •  hr  1 

(  x  1 0 ' 1<,>  ) 

g 

Measured 
Response 
< A  MR  hr  1 

(  x  1 1)  ‘  1  '>  n 

«  K 

^  3  Am 

0.06 

4  2 

20.2 

(0.021) 

20  3 

0  99  5 

l^Cs 

0. 145 

9 . 2 

6  63 
(0.023) 

/  86 

0  849 

137Cs 

0  662 

9.2 

3.86 

(0.024) 

4  03 

0.958 

60Co 

1 .25 

9.2 

3.89 

(0.023) 

4.00 

0.973 

2^3-Am 

0.06 

21 .  7 

20.3 

(0.021) 

21.9 

0.928 

0.145 

21.7 

14.1 

(0.020) 

15.4 

0.916 

137Cs 

0.662 

21.7 

8.57 

(0.022) 

8,78 

0.976 

60Co 

1.25 

21.7 

8.85 

(0.023) 

9.46 

0.936 

24lAm 

0.06 

36.8 

4.05 

(0.022) 

4.92 

0.823 

l^Cs 

0.145 

36.8 

18.1 

(0.021) 

19.3 

0.938 

137Cs 

0.662 

36.8 

12.3 

(0.022) 

12.8 

0.961 

60Co 

1.25 

36.8 

13.1 

(0.023) 

13.2 

0.992 

Experimental  uncertainty 
^Co  source 


±  4%  for 

±  5%  for  34l£s  and  337 c  s  sources 
±  7%  for  241^  source 


Fractional  standard  deviation  (o(x)/x). 


7  7 


k's1'. 


hot  two  cases  which  Involved  low  energy  photon  simtit-s  No  * •  x p  1 . t 
n«»f  ion  could  be  found  as  to  whv  these  two  cases  do  not  agin  ,is  wt  !  1 
is  the  othet  cases  The  discrepancy  for  these  two  cases  could  he 
partially  due  to  an  error  in  t  lie  experiment  (within  the  e  >:pe  i  i  inent  a  1 
procedure  itself  or  in  the  recording  of  the  data)  Furthermore,  the 
consistent  underestimation  of  the  experimental  response  bv  the  MICAP 
calculations  is  not  understood  at  this  time. 

The  comparisons  of  the  MICAP  and  experimental  results  for  the 
AFRRI  50  cm^  ionization  chambers  are  presented  in  Table  4  -  1 5  and  tin- 
comparisons  of  the  MICAP  and  experimental  results  for  the  AFRRI  0 . 5 
cm^  and  0.0b  cm^  ionization  chambers  are  presented  in  Table  4-16. 

The  comparisons  of  calculated  and  experimental  results  for  the  50 
cm3  ionization  chambers  show  excellent  agreement  in  Table  4-15.  The 
calculations  analyzed  500,000  source  particles  incident  on  the 
detector.  The  calculated- to-experimental  ratios  (C/E)  are  well 
within  10%  for  all  cases  analyzed.  The  results  also  show  a  sys¬ 
tematic  overestimation  for  the  C/CO  2  chamber  and  underestimation  for 
the  TE/TE  chamber.  As  in  the  HASL  study,  the  exact  cause  is  not 
known.  A  sensitivity  study  was  performed  however,  on  the  TE/TE  and 
C/CO  2  chambers  with  0.72  cm  wall  thickness.  Both  chambers  were 
evaluated  for  a  ±10%  change  in  gas  density  and  the  C/CO2  chamber  was 
evaluated  for  a  ±5%  change  in  graphite  density.  The  effect  of  gas 
density  was  significant  in  the  sensitivity  calculations.  For  the 
TE/TE  chamber,  the  C/E  ratio  was  0.902  for  a  -10%  change  in  gas  den¬ 
sity  and  1.08  for  a  +10%  change.  Likewise,  for  the  C/CO2  chamber, 
the  C/E  ratio  was  0.877  for  a  -10%  change  in  gas  density  and  1.09 
for  a  +10%  change.  Therefore,  an  error  in  the  gas  density 


Table  4-15.  Comparison  of  Measured  and  Calculated  Detector 
Responses  for  the  Photon  Calibration  Experiments  of  the 
Armed  Forces  Radiobiology  Research  Institute  (AFRRI) 

50  cm^  Ionization  Chambers 


Photon 

Source 

Source 

Energy 

(MeV) 

Detector 

Type 

Wall 

Thickness 

(cm) 

Calculated 

Response 

(pA) 

Measured3 

Response 

(pA) 

C/E 

60Co 

1.25 

C/C02 

0.72 

179  h 

(0.017) 

173 

1.04 

On 

O 

o 

o 

1.25 

C/C02 

1.10 

172 

(0.021) 

171 

1.01 

60Co 

1.25 

C/C02 

1.38 

182 

(0.021) 

170 

1.07 

60Co 

1.25 

C/C02 

1.66 

174 

(0.023) 

168 

1.04 

60Co 

1.25 

TE/TE 

0.72 

123 

(0.017) 

127 

0.969 

60Co 

1.25 

TE/TE 

1.02 

121 

(0.018) 

125 

0.968 

60Co 

1.25 

TE/TE 

1.34 

117 

(0.023) 

124 

0.944 

60Co 

1.25 

TE/TE 

1.67 

120 

(0.023) 

123 

0.976 

Experimental  uncertainty  ±  1% . 


Table  4-16.  Comparison  of  Measured  and  Calculated  Detector 
Responses  for  the  Photon  Calibration  Experiments  of  the 
Armed  Forces  Radiobiology  Research  Institute  (AFRRI) 

0.5  cm^  and  0.05  cm^  Ionization  Chambers 
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calculations  will  result  in  comparable  errors  in  the  calculated 
results.  The  graphite  wall  density  effects  were  found  to  be  less 
significant.  The  sensitivity  study  yielded  a  C/E  ratio  of  1.03  for 
a  -5%  change  in  graphite  density,  and  0.99  for  a  +5%  change.  There¬ 
fore,  a  slight  error  in  graphite  density  will  cause  a  lessor  error 
in  the  calculated  response.  For  the  problems  analyzed,  this  error 
was  within  the  statistical  error  of  the  Monte  Carlo  calculations. 

The  results  of  the  0.5  cm^  and  0.05  cm^  calculations  vary  when 
compared  to  experiment.  As  seen  in  Table  4-16,  the  C/E  ratios  range 
from  0.677  to  1.02  for  the  different  cases  analyzed.  Each  chamber 
was  analyzed  with  900,000  source  particles  incident  on  the  chamber 
wall  material.  Overall,  the  results  indicate  several  trends  in  the 
calculational  ability  of  MICAP  with  respect  to  ionization  chamber 
size.  First,  MICAP  appears  to  have  difficulty  reproducing  the 
smaller  chamber  results.  This  is  evidenced  by  the  0.05  cm^  TE/TE 
C/E  ratio  of  0.677  and  the  0.5  cm^  TE/TE  chamber  (with  no  build-up 
cap)  C/E  ratio  of  0.705.  Variations  in  wall  thickness  (nominally 
0.1  cm)  could  partly  be  responsible  for  these  discrepancies;  how¬ 
ever,  confirmation  would  require  detailed  sensitivity  calculations. 

A  second  trend  is  the  larger  C/E  ratios  for  the  smaller  measured 
signals.  This  could  possibly  be  due  to  errors  in  the  experimental 
results.  It  should  be  noted  that  this  same  trend  is  apparent  in  the 
experimental  uncertainty  presented  at  the  bottom  of  Table  4-16. 

One  final  result  which  contributes  to  the  explanation  of  the 
discrepancy  between  calculated  and  experimental  results  is  seen  in 
the  6®Co  pUise  height  distribution  for  the  50  cm^  TE/TE  ionization 
chamber  shown  in  Figure  4-19.  This  curve  is  representative  of  the 
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Figure  4-19.  Pulse  Height  Distribution  for  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  TE/TE  50  cm^ 

Ionization  Chamber  due  to  a  Mono-Directional  Co  Photon  Source. 


pulse  height  distributions  for  the  ionization  chambers  analyzed  and 


shows  a  broad  distribution  as  a  function  of  energy.  Furthermore, 


the  pulse  height  contains  a  long  high-energy  tail  which  affects  the 


calculation  of  the  average  energy  deposited  per  source  particle. 


Although  the  fractional  standard  deviation  for  the  Monte  Carlo 


result  in  Table  4-15  was  approximately  ±2%,  the  tail  of  the  pulse 


height  distribution  could  be  undersarapled  and  therefore  lead  to  an 


inaccurate  calculated  value. 


The  results  of  the  photon  calibration  experiment  comparisons 


show  the  applicability  of  MICAP  to  ionization  chamber  response  cal¬ 


culations  even  though  the  AFRRI  small  chamber  results  were  incon¬ 


clusive.  The  ultimate  and  final  evaluation  of  MICAP  is  made  with 


the  mixed-field  results  presented  in  the  next  section. 
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4.5  COMPARISON  WITH  MIXED  FIELD  EXPERIMENTS 


The  results  presented  in  the  preceding  section  demons:  r<»t  *•«}  ••• 

ability  of  MICAP  to  calculate  ionization  chamber  response  f unc t  io: • 

for  mono-energetic  photon  calibration  experiments  Although  this 

verification  is  necessary,  the  real  benefit  of  MICAP  is  for  mixed 

field  experiments.  To  demonstrate  the  ability  of  using  MICAP  to 

calculate  ionization  chamber  response  functions  for  mixed  neutron 

and  photon  radiation  environments,  three  radiation  experiments  were 

analyzed. ^2-43  The  fprst  experiment  analyzed  was  performed  at  the 

NBS  252c:f  facility.  This  experiment  involved  exposing  the  50  cm^ 

and  0.5  cm^  ionization  chambers  individually,  in  free  air,  at  a  nom- 

252 

inal  distance  of  30  cm  from  the  center  of  the  Cf  source.  The 
source  and  detector  were  positioned  such  that  room  effects  were 
negligible.  The  second  and  third  experiments  were  performed  at  the 
AFRRI  TRIGA  reactor  described  in  Ref.  44.  These  experiments  were 
performed  in  exposure  room  1  (ER1)  of  the  reactor  facility  for  two 
of  the  most  frequently  used  reactor  configurations: 

1.  Bare  room  with  no  shielding  placed  between  the  experiment 
and  the  reactor  core  (ER1  Free  Field) ,  and 

2.  15  cm  lead  shield  placed  in  front  of  the  reactor  tank  wall 
to  attenuate  photons  (ERl  15  cm  Pb) . 

The  AFRRI  experiments  exposed  a  pair  of  ionization  chambers  about  10 
cm  apart  and  approximately  100  cm  from  the  reactor  core  center.  The 
experimental  setup  for  the  bare  room  (ERl  Free  Field)  is  shown  in 
Figure  4-20  and  for  the  15  cm  lead  shield  (ERl  15  cm  Pb)  in  Figure 
4-21.  The  two  ionization  chambers  exposed  were  the  0.5  cm^  TE/TE 
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Figure  4-20.  Schematic  Top  View  of  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  Experimental  Room  One  Fre 
Field  (ER1  FF)  Experiment. 
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Figure  4-21.  Schematic  Top  View  of  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  Experimental  Room  One  15 
cm  Lead  Shield  (ER1  15  cm  Pb)  Experiment. 


ionization  chamber  and  the  0.5  cra^  Mg/Ar  ionization  chamber.  As  in 
the  NBS  252Cf  experiment,  the  room  effects  were  negligible. 

The  calculations  modeled  the  ionization  chambers  exactly  using 
combinatorial  geometry.  The  material  properties  and  geometry 
description  for  the  AFRRI  chambers  are  given  in  the  last  section  in 
Table  4-13  and  Figures  4-16  and  4-17.  The  primary  neutron  and  pho¬ 
ton  spectra  for  the  AFRRI  exposure  room  experiments  (at  the  detector 
position)  are  given  in  Tables  4-17  and  4-18  respectively.  These 
spectra  were  obtained  from  Ref.  45.  The  NBS  2^2Cf  spectra  (also 
presented  in  Tables  4-17  and  4-18)  were  obtained  from  Refs.  42  and 
46.  The  neutron  source  spectra  (Table  4-17)  was  input  to  the  NEU¬ 
TRON  program  and  the  recoil  heavy  ion  and  secondary  photon  produc¬ 
tion  sources  were  generated  for  subsequent  transport  in  program 
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Table  4-17.  Primary  Neutron  Source  Spectra  Used  in  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  Mixed  Neutron  and  Phouon 
Radiation  Field  Ionization  Chamber  Calculations 


Group 

Number 

Neutron 
Energy  (eV) 

ERl* 

Free  Field 

ERl* 

15  cm  Lead 

252Cf* 

Source 

1 

1.96+07 

4.95-06 

2 . 14-06 

1.62-05 

2 

1.69+07 

2.25-05 

9.42-06 

5.88-05 

3 

1.49+07 

2.12-05 

8 . 59-06 

4.94-05 

4 

1.42+07 

1.52-05 

6.00-06 

4.12-05 

5 

1.38+07 

8.80-05 

3.46-05 

1.71-04 

6 

1.28+07 

7.79-05 

3 . 00-05 

1 . 76-04 

7 

1.22+07 

3 . 13-04 

1.19-04 

5.92-04 

8 

1.11+07 

6.68-04 

2.55-04 

1.26-03 

9 

1.00+07 

1.25-03 

4.73-04 

2.19-03 

10 

9.05+06 

2.06-03 

7.66-04 

3.65-03 

11 

8 . 19+06 

3.30-03 

1.23-03 

5.73-03 

12 

7.41+06 

8.19-03 

3.04-03 

1.39-02 

13 

6.38+06 

2.49-02 

8.82-03 

4.24-02 

14 

4.97+06 

6.55-03 

2.02-03 

1.21-02 

15 

4.72+06 

2.21-02 

6.85-03 

4.15-02 

16 

4.07+06 

9 . 74-02 

2 . 70-02 

1.17-0; 

17 

3.01+06 

8.40-02 

3.68-02 

1.02-01 

18 

2.39+06 

8 . 80-03 

1.03-02 

1.57-02 

19 

2.31+06 

8.46-02 

5.40-02 

1.08-01 

20 

1.83+06 

1.17-01 

1 . 60-01 

2.06-01 

21 

1.11+06 

1 . 18-01 

2.14-01 

1.87-01 

22 

5.50+05 

9.60-02 

1.69-01 

1.23-01 

23 

1.58+05 

2.18-02 

3.20-02 

9.66-03 

24 

1.11+05 

2.97-02 

5 . 14-02 

9.03-03 

25 

5.25+04 

2.29-02 

3.14-02 

2.85-03 

26 

2.48+04 

4.74-03 

2.80-03 

2 . 30-04 

27 

2 . 19+04 

2.62-02 

1.73-02 

7 . 56-04 

28 

1.03+04 

3.55-02 

3.44-02 

2.91-04 

29 

3.35+03 

2.70-02 

1.45-02 

5.13-05 

30 

1.23+03 

2.06-02 

1.66-02 

9.87-06 

31 

5.83+02 

5.10-02 

2.86-02 

4.44-06 

32 

1.01+02 

2.81-02 

2.73-02 

2.92-07 

33 

2 . 90+01 

2.05-02 

1 . 34-02 

4.11-08 

34 

1.07+01 

1.31-02 

1.46-02 

1.01-08 

35 

3.06+00 

8.40-03 

8.63-03 

1.41-09 

36 

1.13+00 

6.02-03 

5.03-03 

3.17-10 

37 

4 . 14-01 
1.00-05 

9.74-03 

7.61-03 

9.04-11 

*Flux  spectra  presented  as  normalized  probability  distribution 
functions . 
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Table  4-18.  Primary  Photon  Source  Spectra  Used  in  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRI)  Mixed  Neutron  and  Photon 
Radiation  Field  Ionization  Chamber  Calculations 


*Flux  spectra  presented  as  normalized  probability  distribution 
functions . 


HEAVY  and  program  PHOTON,  respectively.  Because  the  detector  effi¬ 
ciency  (with  respect  to  neutron  interactions)  was  so  poor,  a  ficti¬ 
tious  scattering  model  was  incorporated  in  the  NEUTRON  program  to 
facilitate  interactions.  This  eliminated  the  capability  of  produc¬ 
ing  a  pulse  height  distribution,  and  therefore  one  analog  Monte 
Carlo  case  was  also  performed.  With  the  fictitious  scattering 
model,  approximately  200,000  to  500,000  source  particles  were 
analyzed  to  obtain  a  fractional  standard  deviation  less  than  ±10^. 
These  cases  used  approximately  15  to  35  minutes  of  cpu  time  on  an 
IBM  3033.  The  analog  case  modeled  3,000,000  source  particles  to 
obtain  the  same  fractional  standard  deviation  and  used  approximately 
3  hours  and  20  minutes  cpu  time.  As  in  the  case  with  the  photon 
calibration  experiments,  the  primary  photon  source  spectra  (Table 
4-18)  were  modeled  in  program  PHOTPREP  to  generate  an  input  source 
tape  for  program  PHOTON.  The  PHOTON  source  tape  modeled  900,000 
photons  incident  on  the  ionization  chamber  wall.  The  results  from 
program  HEAVY  and  program  PHOTON  were  analyzed  to  obtain  the  total 
energy  deposited  in  the  cavity  by  particle  type  i.e.,  protons,  alpha 
particles,  primary  photons,  etc.  These  results  were  then  normalized 
to  the  incident  source  strength  to  obtain  the  neutron,  photon  and 
total  response  of  the  ionization  chamber.  The  cross  sections  and 
material  data  for  the  transport  programs  were  generated  using  the 
neutron  and  photon  cross  section  processor  modules  MICRO,  MACRO,  and 
PECSP.  Recombination  effects  were  assumed  negligible,  and  therefore 
were  not  considered. 

The  comparisons  of  the  calculated  and  experimental  total  ioni¬ 
zation  chamber  responses  to  the  mixed  field  experiments  are 


presented  In  Table  4-19.  The  calculated- to -experimental  ratios 
(C/E)  show  excellent  agreement  in  all  the  cases  except  the  0.5  cm 
Mg/Ar  chamber  in  the  ER1  15  cm  Pb  field.  Like  the  mono-energetic 
photon  results,  the  largest  discrepancy  (excluding  the  case  men¬ 
tioned  above)  occurs  for  the  smallest  measured  response.  A  portion 
of  this  discrepancy  could  therefore  be  due  to  experimental  error. 

The  Monte  Carlo  results  have  a  fractional  standard  deviation  less 
than  ±10%  for  all  calculations.  This  fractional  standard  deviation 
is  for  the  detector  response  per  incident  source  particle  i.e., 

Amps/ (n-cm^- sec -kW)  or  Amps/(7-cm^-sec-kW) .  To  compare  with  the 
measured  response,  an  absolute  normalization  was  obtained  from  the 
spectra  presented  in  Ref.  45  for  the  two  AFRRI  experiments.  While 
the  spectral  shapes  in  Ref.  45  are  probably  accurate  to  within  ±10%, 
the  absolute  normalization  is  known  to  within  only  ±15-20%,  espe¬ 
cially  for  the  AFRRI  ER1  15  cm  Pb  cases.  Therefore,  the  MICAP  cal¬ 
culated  results  in  Table  4-19  are  within  the  statistical  error  of 
the  information  used. 

The  AFRRI  ER1  15  cm  Pb  calculation  using  the  0.5  cm^  Mg/Ar  ion¬ 
ization  chamber  displayed  an  interesting  result  in  the  neutron 
response  analysis.  In  viewing  the  geometry  description  in  Figure 
4-17,  a  thin  polyethylene  strip  (approximately  0.1-cm-high)  is  shown 
in  contact  with  the  gas  region.  The  polyethylene  acts  as  an  insula¬ 
tor  between  the  outside  wall  and  the  central  electrode.  The 
analysis  of  the  results  from  program  HEAVY  shows  recoil  argon  ions 
from  neutron  interactions  in  the  gas,  and  proton  recoils  from  neu¬ 
tron  interactions  with  hydrogen  in  the  polyethylene  strip,  as  the 
only  two  sources  of  energy  deposition  in  the  gas.  Furthermore. 


Table  4-19.  Comparison  of  Measured  and  Calculated  Detector 
Responses  for  the  Mixed  Neutron  and  Photon  Radiation 
Field  Experiments  of  the  Armed  Forces  Radiobiology 
Research  Institute  (AFRRI)  Ionization  Chambers 


Detector  Calculated  Measured 

Radiation  Volume  Detector  Response  Response 


Field 

(cmJ) 

Type 

(pA) 

(pA) 

C/E 

252cf 

50 

C/C02 

6.51 

(0.023) 

7.21 

0.903 

252Cf 

50 

TE/TE 

12.2 

(0.062) 

12.4 

0.984 

252Cf 

0.5 

Mg/Ar 

0.061 

(0.019) 

0.078 

0.782 

252Cf 

0.5 

TE/TE 

0.105 

(0.012) 

0.130 

0.808 

ER1  FF 

0.5 

Mg/Ar 

7.10 

(0.037) 

6.42 

1.11 

ER1  FF 

0.5 

TE/TE 

8.21 

(0.021) 

7.70 

1.07 

ER1  15cm 

Pb 

0.5 

Mg/Ar 

4.42 

(0.072) 

7.16 

0.617 

ER1  15cm 

Pb 

0.5 

TE/TE 

8.66 

8.75 

0.990 

(0.007) 


Experimental  uncertainty  ±  5%. 
^Fractional  standard  deviation  (<7(x)/x)  . 


approximately  80%  of  the  total  signal  (due  to  neutrons)  originates 
from  the  proton  recoil  out  of  the  polyethylene  strip.  This  result 
is  significant  in  that  undersampling  has  probably  occurred  in  the 
MICAP  results  for  the  polyethylene  strip  in  the  0.5  cm^  Mg/Ar 
chamber.  Because  the  photon  component  of  the  total  signal  is  dom¬ 
inant  in  the  ^^Cf  ancj  t^e  afrri  ER1  FF  results,  the  undersam¬ 
pling  would  not  cause  a  large  discrepancy  in  a  comparison  of  calcu¬ 
lated  and  experimental  total  responses.  The  AFRRI  ER1  15  cm  Pb 
field  is  predominantly  a  neutron  field  however,  and  the  0.5  cm^ 

Mg/Ar  ionization  chamber  is  designed  to  be  highly  neutron  insensi¬ 
tive.  The  total  response  of  the  0.5  cm^  Mg/Ar  ionization  chamber 
remains  dominated  by  the  neutron  component,  and  therefore  an  under¬ 
sampling  (and  subsequent  underestimation)  of  the  neutron  component 
would  cause  a  significantly  lower  calculated  value.  The  results 
appear  to  justify  this  scenario  with  a  C/E  of  0.617.  A  possible 
solution  which  surfaces  from  the  above  argument  would  be  to  replace 
the  polyethylene  with  a  non-hydrogeneous  insulating  material.  This 
would  mitigate  the  proton  recoil  component  of  the  neutron  signal  and 
make  the  0.5  cm^  Mg/Ar  chamber  truly  neutron  insensitive. 

The  experimental  results  use  the  "two  dosimeter"  method 
described  in  Ref.  47  to  determine  the  neutron  and  photon  contribu¬ 
tion  to  the  total  signal.  Briefly,  the  method  requires  two  dosime¬ 
ters  with  different  neutron  sensitivities.  One  dosimeter,  such  as 
the  TE/TE  ionization  chamber,  will  have  approximately  the  same  sen¬ 
sitivity  to  neutrons  and  photons.  The  second  dosimeter,  such  as  the 
Mg/Ar  ionization  chamber  will  have  a  reduced  neutron  sensitivity 
relative  to  photons.  The  neutron  and  photon  components  at  a  point 


in  the  radiation  field  can  be  computed  from  the  response  of  the  two 
ionization  chambers  through  the  use  of  two  simultaneous  equations 
and  the  detector  responses  to  the  photon  calibration  experiments. 
MICAP  on  the  other  hand,  requires  only  the  analysis  of  the  ioniza¬ 
tion  chamber  with  approximately  equal  neutron  and  photon  sensitivi¬ 
ties  to  ascertain  the  relative  components  of  the  total  response.  A 
comparison  of  these  results  in  Table  4-20  show  excellent  agreement 
thus  reaffirming  the  ability  of  MICAP  to  calculate  mixed  field  radi¬ 
ation  experiments. 

The  final  results  presented  for  the  mixed  neutron  and  photon 
radiation  fields  are  the  pulse  height  distributions  for  the  50  cm^ 
TE/TE  ionization  chamber  in  the  ^52cf  source  field.  As  stated  ear¬ 
lier,  only  one  analog  Monte  Carlo  case  was  performed  and  therefore 
only  one  set  of  pulse  heights  are  possible.  The  pulse  height  dis¬ 
tribution  due  to  the  ^52Cf  photon  source  presented  in  Figure  4-22  is 
very  similar  to  the  ^Co  pUise  height  presented  in  Figure  4-19. 
Therefore,  the  arguments  presented  with  respect  to  the  ^®Co  pulse 
height  distribution  are  equally  applicable  here.  The  pulse  height 
distribution  for  the  252^f  neutron  source  presented  in  Figure  4-23 


depicts  a  continuously  decreasing  curve  as  a  function  of  increasing 
energy.  While  there  are  significantly  more  contributions  in  the 
lower  energy  bins,  these  particles  are  outweighed  by  much  larger 
energy  depositions  of  the  fewer  particles  in  the  higher  energy  bins. 
For  the  pulse  height  distribution  shown  in  Figure  4-23,  the  average 
energy  deposited  is  approximately  3.0  x  10'^  MeV/neutron.  As  in  the 
case  with  the  photon  pulse  height  distributions,  undersampling  of 
the  long  high  energy  tail  in  the  neutron  pulse  height  distribution 
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could  lead  to  an  underestimation  of  the  average  energy  deposited  and 
therefore  may  lead  to  an  inaccurate  calculated  value. 


Table  4-20.  Comparison  of  Measured  and  Calculated  Detector  Response 
Components  for  the  Mixed  Neutron  and  Photon  Radiation 
Field  Experiments  of  the  Armed  Forces  Radiobiology 
Research  Institute  (AFRRI)  Ionization  Chambers 


Radiation 

Field 

Dosimeter 

Pair 

Calculated  Response 
Photon/Total 

Measured  Response 
Photon/Total 

252Cf 

te/te-c/co2 

0.31 

(0.064) 

0.33 

252Cf 

TE/TE-Mg/Ar 

0.26 

(0.025) 

0.31 

ER1  FF 

TE/TE-Mg/Ar 

0.63 

(0.012) 

0.66 

ERl  15cm  Pb 

TE/TE-Mg/Ar 

0.10 

(0.043) 

0.09 

Experimental  uncertainty  ±  5% . 
^Fractional  standard  deviation  (cr(x)/x) . 
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Figure  4-22.  Pulse  Height  Distribution  for  the  Armed  Forces 
Radiobiology  Research  Institute  (AFRRX)  TE/TE  50  cm^  Ionization 

252 

Chamber  due  to  a  Mono-Directional  Cf  Neutron  Source. 


OUL  DUC  W  l»M 


ENERGY  (MeV) 

Figure  4-23.  Pulse  Height  Distribution  for  the  Armed  Forces 
Radloblology  Research  Institute  (AFRRI)  TE/TE  50  cm^  Ionization 

252 

Chamber  due  to  a  Mono-Directional  Cf  Photon  Source. 
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CHAPTER  V 


CONCLUSIONS 


This  research  project  has  consisted  of  the  development  and 


verification  of  a  new  Monte  Carlo  code  system  for  calculating 


ionization  chamber  responses  to  mixed  neutron  and  photon  radiation 


environments.  This  new  code  system,  MICAP  -  Monte  Carlo  Ionization 


Chamber  Analysis  Package,  represents  a  unique  capability  in  the  area 


of  ionization  chamber  response  analysis  for  mixed  field  dosimetry. 


A  five  pa^t  verification  program  was  used  to  examine  all  aspects  of 


the  radiation  processes  applicable  to  ionization  chamber  analysis. 


The  verification  program  involved  comparisons  with  results  obtained 


from  existing  codes  as  well  as  comparisons  to  experimental  data. 


The  comparisons  with  existing  codes  verified  the  individual 


aspects  of  MICAP  with  respect  to  ionization  chamber  analysis.  More 


specifically,  iron  transmission  comparisons  with  MORSE  yielded 


excellent  agreement  (within  the  Monte  Carlo  statistical  error)  in 


the  neutron  leakage  spectra,  secondary  photon  production  spectra, 


and  the  photon  leakage  spectra  for  various  neutron  sources  and  iron 


slab  thicknesses.  Secondly,  comparisons  with  MACK-IV  and  RECOIL 


data  verified  the  nuclear  models  used  in  MICAP  to  describe  the 


individual  nuclear  reactions  with  an  agreement  to  within  10%  for 


average  energy  depositions  and  kerma  factors.  Final lv.  the 


comparison  to  05S ,  a  specialized  Monte  Carlo  detector  analysis  cod« 


which  has  been  used  extensively,  showed  excellent  agreement  in  tin 


pulse  height  distributions  for  a  plastic  scintillator,  and 


therefore,  verified  the  charged  particle  transport  \ 


tract-  skis 


P 


I 


I 


The  comparisons  with  experimental  data  first  evaluated  MICAP 
with  respect  to  photon  calibration  experiments.  The  results  from 
these  comparisons  yielded  agreement  generally  to  within  10%  for  the 
ionization  chamber  responses  to  the  mono-energetic  photon  sources. 
Finally,  the  comparisons  with  experiments  performed  in  the  mixed 
neutron  and  photon  radiation  environments  at  NBS  and  AFRRI  verified 
the  applicability  of  MICAP  to  mixed  field  ionization  chamber 
response  analysis.  The  calculated  results  agreed  with  experimental 
results  to  within  the  statistical  error  of  the  Monte  Carlo  and  the 
normalization  of  the  experimental  data. 

The  final  product  (MICAP)  represents  a  valuable  tool  to 
radiation  dosimetry.  The  modular  code  system  provides  the  most 
rigorous  treatment  to  date  of  the  particle  transport  processes 
applicable  to  ionization  chamber  response  analysis.  The  general 
implementation  procedures  incorporated  into  MICAP  allow  the  system 
to  be  useful  for  many  applications  in  the  field  of  radiation 
dosimetry.  It  is  recommended  that  the  code  system  be  tested  further 
by  additional  comparisons  with  mixed  field  radiation  experiments  to 
resolve  the  discrepancies  noted  in  this  document.  Additional 
recommendations  for  future  work  are  presented  in  the  next  chapter. 
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CHAPTER  VI 


RECOMMENDATIONS  FOR  FUTURE  WORK 

This  study  has  demonstrated  the  applicability  of  using  a  gen¬ 
eral  Monte  Carlo  code  system  for  calculating  an  ionization  chamber 
response  to  a  mixed  neutron  and  photon  radiation  environment.  The 
FORTRAN  code  system  developed  and  tested  in  this  study  is  only  a 
prototype,  not  a  production  code.  The  current  version  of  MICAP 
could  be  improved  along  several  lines  within  the  present  scope  of 
ionization  chamber  response  calculations.  In  particular,  the  fol¬ 
lowing  modifications  and/or  investigations  of  MICAP  are  suggested: 

1.  Improving  the  programming  through  the  use  of  more  effi¬ 
cient  algorithms,  ASSEMBLER  language  coding  of  the  most 
cpu- intensive  subroutines,  and  modifying  program  input  to 
establish  a  more  user  friendly  code  system. 

2.  Incorporating  more  sophisticated  and  detailed  charged  par 
tide  evaporation  models  for  the  neutron  induced  charged 
particle  production  reactions  i.e.,  (n,p),  (n,a),  etc. 

The  more  sophisticated  models  should  yield  a  better 
representation  of  the  energy  of  the  charged  particles  pro 
duced  via  neutron  interactions. 

3.  Performing  more  comparisons  with  mixed  neutron  and  photon 
radiation  experiments  -  especially  experiments  involving 
high  energy  neutron  sources  to  test  the  accuracy  of  the 


charged  particle  evaporation  models  and  further  substan¬ 
tiate  the  applicability  of  MICAP  to  ionization  chamber 
response  analysis. 

4.  Study  the  incorporation  of  a  multiple  scattering  capabil¬ 
ity  in  the  low  energy  charged  particle  transport  program 
using  theoretical  and  empirical  formulas.  The  current 
work  contains  only  continuous  energy  loss  models.  Allow¬ 
ing  multiple  scattering  in  the  low  energy  charged  particle 
transport  will  yield  a  more  accurate  representation  of  the 
physics . 

5.  Incorporating  the  fictitious  scattering  model  used  in  the 
present  work  as  an  input  option  to  facilitate  calculations 
involving  ionization  chambers  which  exhibit  low  neutron 
efficiency. 

6.  Completing  the  programming  for  the  Monte  Carlo  source 
biasing  to  improve  the  calculation  statistics. 

Finally,  MICAP  could  be  further  generalized  to  operate  as  a  point 
Monte  Carlo  code  system  applicable  to  reactor  core  and  reactor 
shielding  calculations.  This  could  be  accomplished  through  program 
switches  which  engage  (or  dissengage)  the  portions  of  the  code  sys¬ 
tem  applicable  to  a  particular  problem.  Such  a  code  system  would 
truly  represent  a  general  Monte  Carlo  program  applicable  to  all 
aspects  of  radiation  transport/dosimetry  analysis. 


98 


LIST  OF  REFERENCES 


1.  H.  W.  FuLbright,  "Ionization  Chambers,"  Nuclear  Instruments  and 
Methods,  Vol  162,  Nos.  1-3,  Part  1,  pp.  21-28,  1979. 

2.  T.  E.  Burlin,  "Cavity  Chamber  Theory,"  Radiation  Dosimetry , 

Vol.  1,  F.  H.  Attix,  W.  C.  Roesch,  and  E.  Tochlin,  Eds., 
Academic  Press,  New  York,  1968. 

3.  Radiation  Dosimetry:  X  Rays  and  Gamma  Rays  with  Maximum  Photon 
Energies  Between  0.6  and  50  MeV ,  ICRU  Report  14,  International 
Committee  on  Radiation  Units  and  Measurements,  Washington, 

D.C.,  1969. 

4.  Radiation  Dosimetry:  X  Rays  Generated  at  Potentials  of  5  to 
150  keV ,  ICRU  Report  17,  International  Committee  on  Radiation 
Units  and  Measurements,  Washington,  D.C.,  1970. 

5.  Measurement  of  Absorbed  Dose  in  a  Phantom  Irradiated  by  a  Sin¬ 
gle  Beam  of  X  or  Gamma  Rays,  ICRU  Report  23,  International  Com¬ 
mittee  on  Radiation  Units  and  Measurements,  Washington,  D.C., 
1973. 

6.  H.  Schraube,  "Concluding  Discussion  on  Ion  Chambers  for  Neutron 
Dosimetry,"  Ion  Chambers  for  Neutron  Dosimetry ,  J.  J.  Broerse, 
Ed.,  EUR  6782,  Commission  of  the  European  Community,  Luxem¬ 
bourg,  pp.  345-350,  1980. 

7.  An  International  Neutron  Dosimetry  Intercomparison  (INDI) ,  ICRU 
Report  27,  International  Committee  on  Radiation  Units  and  Meas¬ 
urements,  Washington,  D.C.,  1978. 

8.  J.  J.  Broerse,  G.  Burger,  and  M.  Coppola,  A  European  Neutron 
Dosimetry  Intercomparison  Project  (ENDIP) ,  EUR  6004,  Commission 
of  the  European  Community,  Luxembourg,  1978. 

9.  H.  Bichsel,  J.  Eenmaa,  K.  Weaver,  and  P.  Wooten,  "Dosimetry  of 
Mixed  Radiation  Fields,"  proceedings  of  a  Symposium  on  Biologi¬ 
cal  Dosimetry,  IAEA-SM-193/78 ,  pp.  3-17,  Int.  Atomic  Energy 
Agency,  Vienna,  1975. 

10.  Neutron  Dosimetry  for  Biology  and  Medicine,  ICRU  Report  26, 
International  Committee  on  Radiation  Units  and  Measurements, 
Washington,  D.C.,  1977. 

11.  H.  Ing  et  al..  Compendium  of  Neutron  Spectra  in  Criticality 
Accident  Dosimetry ,  IAEA- 180,  Int.  Atomic  Energy  Agency, 

Vienna,  1978. 


100 


■ 


12.  J.  E.  Turner  et  al.,  User's  Manual  for  PHOEL-2,  A  Monte  Carlo 
Computer  Code  for  Calculating  Energies  of  Photoelectrons  and 
Compton  Electrons  in  Water,  ORNL/TM-6954,  Union  Carbide  Corp. 
Nuclear  Div. ,  Oak  Ridge  Natl.  Lab.,  1979. 

13.  J.  E.  Turner  et  al.,  "Calculated  Electron  Slowing  Down  Spectra 
for  Liquid  Water  Irradiated  by  X  and  Gamma  Rays-Implications 
for  Photon  RBE,"  Sixth  Symposium  on  Microdosimetry ,  J.  Booz  and 
H.  G.  Ebert,  Eds.,  Harwood  Academic  Publishers,  London,  1978. 

14.  S.  Pszona  and  M.  Makarewicz ,  "Effect  of  Cavity  Size  on  the  Sen¬ 
sitivity  of  a  TE-Walled,  TE-Gas- Filled  Ionization  Chamber  for 
Fast  Neutrons,"  Physics  in  Medicine  and  Biology,  Vol.  27,  No. 

8,  pp.  1015-1022,  1982. 

15.  M.  W.  Waddell,  Jr.,  User's  Guide  to  PXMORSE ,  ORNL/CSD/TM-129 , 
Union  Carbide  Corp.  Nuclear  Di/.,  Oak  Ridge  Natl.  Lab.,  1980. 

16.  M.  W.  Waddell,  Jr.,  "A  Continuous  Energy  Monte  Carlo  Code," 

M.S.  thesis,  The  University  of  Tennessee,  Knoxville,  1980. 


17.  M.  B.  Emmett,  The  MORSE  Monte  Carlo  Radiation  Transport  Code 
System,  ORNL-4972,  Union  Carbide  Corp.  Nuclear  Div.,  Oak  Ridge 
Natl.  Lab.,  1975. 

18.  T.  A.  Gabriel,  J.  D.  Amburgey,  and  B.  L.  Bishop,  CALOR-A  Monte 
Carlo  Program  Package  for  the  Design  and  Analysis  of  Calorime¬ 
ter  Systems,  ORNL/TM-5619 ,  Union  Carbide  Corp.  Nuclear  Div. , 
Oak  Ridge  Natl.  Lab.,  1977. 

19.  R.  L.  Ford  and  W.  R.  Nelson,  The  EGS  Code  System :  Computer 
Program  for  the  Monte  Carlo  Simulation  of  Electromagnetic  Cas¬ 
cade  Showers,  SLAC-210,  Stanford  Linear  Accelerator  Center, 
1978. 

20.  Average  Energy  Required  to  Produce  an  Ion  Pair,  ICRU  Report  31, 
International  Committee  on  Radiation  Units  and  Measurements, 
Washington,  D.C.,  1979. 


L? 

K 

V. 

F 

y 

»• 


J.  B.  Birks,  The  Theory  and  Practice  of  Scintillation  Counting , 
The  Macmillian  Company,  New  York,  1964. 


Los  Alamos  Radiation  Transport  Group  X-6,  MCNP-A  Generalized 
Monte  Carlo  Code  for  Neutron  and  Photon  Transport ,  LA- 7396 -M, 
Los  Alamos  Natl.  Lab.,  revised  April  1981. 


J.  C.  Nimal  et  al.,  Programme  de  Monte  Carlo  Polycinetique  a 
Trois  Dimensions,  TRIPOLI-2,  C.E.A.,  Saclay,  1980. 


S.  N.  Cramer  et  al.,  The  Use  of  Kernels  in  Studying  Neutron 
Transport  Problems ,  ORNL-TM- 2508 ,  Union  Carbide  Corp.  Nuclear 
Div.,  Oak  Ridge  Natl.  Lab.,  1969. 


fll 


bl 


D.  E.  Bartine  et  al . ,  Low-Energy  Electron  Transport  by  the 
Method  of  Discrete  Ordinates ,  ORNL-TM-3438,  Union  Carbide  Corp. 
Nuclear  Div. ,  Oak  Ridge  NaCl.  Lab.,  1971. 

R.  Kinsey,  ENDF-102  Data  Formats  and  Procedures  for  the 
Evaluated  Nuclear  Data  File,  ENDF/B-V,  BNL-NCS- 50496, 

Brookhaven  Natl.  Lab.,  1979. 

L.  L.  Carter  and  E.  D.  Cashwell,  Particle  Transport  Simulation 
with  the  Monte  Carlo  Method,  TID-26607,  USERDA  Technical  Infor¬ 
mation  Center,  1975. 

D.  G.  Doran  and  N.  J.  Graves,  Displacement  Cross  Sections  and 
PKA  Spectra:  Tables  and  Applications ,  HEDL-TME  76-70,  Hanford 
Natl.  Lab.,  1976. 

L.  Dresner,  EVAP-A  FORTRAN  Program  for  Calculating  the  Evapora¬ 
tion  of  Various  Particles  from  Excited  Compound  Nuclei,  ORNL- 
TM-196,  Union  Carbide  Corp.  Nuclear  Div.,  Oak  Ridge  Natl.  Lab., 
1961. 

P.  N.  Stevens  and  H.  C.  Claiborne,  Weapons  Radiation  Shielding 
Handbook ,  DASA-1982-5,  Defense  Atomic  Support  Agency,  Washing¬ 
ton,  D.C. ,  1970. 

N.  M.  Schaeffer,  Reactor  Shielding  for  Nuclear  Engineers ,  TID- 
25951,  USERDA  Technical  Information  Center,  1973. 

T.  W.  Armstrong  and  K.  C.  Chandler,  SPAR  A  FORTRAN  Program  for 
Computing  Stopping  Powers  and  Ranges  for  Muons,  Charged  Pions, 
Protons,  and  Heavy  Ions,  ORNL-4869,  Union  Carbide  Corp.  Nuclear 
Div.,  Oak  Ridge  Natl.  Lab.,  1973. 

J.  0.  Johnson  and  T.  A.  Gabriel,  A  User's  Guide  to  MICAP  -  A 
Monte  Carlo  Ionization  Chamber  Analysis  Package,  ORNL/TM- 10340 , 
Martin  Marietta  Energy  Systems,  Inc.,  Oak  Ridge  Natl.  Lab.,  (to 
be  published  in  1987). 

M.  A.  Abdou,  Y.  Gohar,  and  R.  Q.  Wright,  MACK - IV ,  A  New  Version 
of  MACK :  A  Program  to  Calculate  Nuclear  Response  Functions  for 
Data  in  ENDF/B  Format,  ANL/FPP/TM-77-5,  Argonne  Natl.  Lab., 
1978. 

T.  A.  Gabriel,  J.  D.  Amburgey,  and  N.  M.  Greene,  Radiation- 
Damage  Calculations :  Primary  Recoil  Spectra,  Displacement 
Rates,  and  Gas -Product ion  Rates,  0RNL/TM-5160,  Union  Carbide 
Corp.  Nuclear  Div.,  Oak  Ridge  Natl.  Lab.,  1976. 

R.  E.  Textor  and  V.  V.  Verbinski,  OSS :  A  Monte  Carlo  Code  for 
Calculating  Pulse  Height  Distributions  Due  to  Monoenergetic 
Neutrons  Incident  on  Organic  Scintillators ,  ORNL-4160,  Union 
Carbide  Corp.  Nuclear  Div.,  Oak  Ridge  Natl.  Lab.,  1968. 


37. 

! 

S.  N.  Cramer  and  R.  W.  Roussln,  Monte  Carlo  Analysis  of  a  1 
Time -Dependent  Neutron  and  Secondary  Gamma  Ray  Integral  Experi-  1 
ment  on  a  Thick  Concrete  and  Steel  Shield,  ORNL/TM- 7972 ,  Union  ■ 
Carbide  Corp.  Nuclear  Div. ,  Oak  Ridge  Natl.  Lab.,  1981.  E 

- 

38. 

D.  E.  Bartine  et  al . ,  37  Neutron,  21  Gamma  Ray  Coupled,  P3 , 

Multigroup  Library  in  ANISN  Format,  DLC- 31/(DPL- 1/FEWG1 ) ,  1975.  J 

39. 

j 

J.  A.  DeCampo,  H.  L.  Beck,  and  P.  D.  Raft,  High  Pressure  Argon  g 
Ionization  Chamber  Systems  for  the  Measurement  of  Environmental  jj 
Radiation  Exposure  Rates,  HASL-260,  Health  and  Safety  Labora-  X 
tory,  New  York,  1972. 

40. 

1  Lt.  M.  A.  Dooley,  Armed  Forces  Radiobiological  Research  a 
Institute,  letter  to  J.  0.  Johnson,  Oak  Ridge  Natl.  Lab.,  Oak  | 
Ridge,  Tennessee,  dated  May  18  1984.  1 

i 

41. 

1  Lt.  M.  A.  Dooley  et  al . ,  "Wall  Attenuation  and  Scatter 

Characteristics  of  Ionization  Chambers  at  the  Armed  Forces 

Radiobiology  Research  Institute,"  draft  report,  Armed  Forces  i 

Radiobiology  Research  Institute,  Bethesda,  Maryland,  1986.  1 

» 

42. 

D.  E.  Eagleson,  Armed  Forces  Radiobiology  Research  Institute  J 
(AFRRI),  Bethesda,  Maryland,  private  communications  to  J.  0.  ij 
Johnson,  Oak  Ridge  Natl.  Lab.,  Oak  Ridge,  Tennessee,  May-July  j 
1986.  jj 

0 

43. 

1  Lt.  M.  A.  Dooley  et  al.,  "NBS-AFRRI  Ionization  Chamber  Inter-  J 
comparison  in  Mixed  Neutron  and  Gamma-Ray  Radiation  Fields,"  3 
TR-86-3,  Armed  Forces  Radiobiology  Research  Institute,  1986.  j 

44. 

J.  A.  Sholtis,  Jr.  and  M.  L.  Moore,  Reactor  Facility ,  Armed  2 
Forces  Radiobiology  Research  Institute ,  AFRRI  TR  81-2,  Armed  1 
Forces  Radiobiology  Research  Institute,  Bethesda,  Maryland,  3 
1981.  jj 

45. 

V.  V.  Verbinski  and  C.  G.  Cassapakis,  Calculation  of  the  Neu-  ^ 
tron  and  Gamma-Ray  Environment  In  and  Around  the  AFRRI  TRIGA  j 
Reactor,  DNA  5793F-2,  Defense  Nuclear  Agency,  Washington,  D.C.,  j 
198l.  j 

46. 

i 

C.  M.  Eisenhauer,  National  Bureau  of  Standards,  Washington,  D.  j 
C.,  private  communication  to  J .  0.  Johnson,  Oak  Ridge  Natl.  a 
Lab.,  Oak  Ridge,  Tennessee,  July  1986.  J 

• 

47. 

l 

American  Association  of  Physicists  in  Medicine,  "Protocol  for 

Neutron  Beam  Dosimetry,"  AAPM  Report  No.  7  of  Task  Group  18, 

Fast  Neutron  Beam  Physics,  Radiation  Therapy  Committee,  1980. 

f 

» 

1 

103  J 

i 

I 

» 

i 

.  _  _  _ _ _  ..... _ ...J 

ORNL/TM- 10196 


INTERNAL  DISTRIBUTION 


1. 

F.  S.  Alsmiller 

32. 

R.  T.  Santoro 

2. 

R.  G.  Alsmiller,  Jr. 

33. 

M.  S.  Smith 

3. 

D.  E.  Bartine 

34. 

J .  S .  Tang 

4. 

B.  L.  Bishop 

35. 

M.  W.  Waddell 

5. 

T .  J .  Burns 

36. 

R.  C.  Ward 

6. 

D.  G.  Cacuci 

37. 

R.  M.  Westfall 

7. 

D.  H.  Cook 

38. 

A.  Zucker 

8. 

S.  N.  Cramer 

39. 

P.  W.  Dickson,  Jr. 

9. 

J.  K.  Dickens 

(Consultant) 

10. 

J.  D.  Drischler 

40. 

G.  H.  Golub 

11. 

M.  B.  Emmett 

(Consultant) 

12. 

W.  W.  Engle,  Jr. 

41. 

R.  Haralick 

13 

-17. 

T.  A.  Gabriel 

(Consultant) 

18. 

D.  T.  Ingersoll 

42. 

D.  Steiner 

19 

-23. 

J.  0.  Johnson 

(Consultant) 

24. 

R.  A.  Lillie 

43. 

Central  Research  Library 

25. 

F.  C.  Maienschein 

44-48. 

EPMD  Reports  Office 

26. 

J.  V.  Pace  III 

49. 

ORNL  Y-12  Technical  Library 

27. 

R.  W.  Peelle 

Document  Reference  Section 

28. 

F.  G.  Perey 

50-51. 

Laboratory  Records 

29. 

L.  M.  Petrie 

52. 

ORNL  Patent  Office 

30. 

R.  W.  Roussin 

31. 

J .  C .  Ryman 

EXTERNAL  DISTRIBUTION 

53.  Dr.  Charles  Aeby,  Air  Force  Weapons  Laboratory,  Kirkland  Air  Force 
Base,  Albuquerque,  NM  87117. 

54.  Office  of  Assistant  Manager  for  Energy  Research  &  Development, 
DOE-ORO,  Oak  Ridge,  TN  37830. 

55.  Argonne  National  Laboratory,  Library  Services  Department,  302-CE125, 
9700  S.  Cass  Avenue,  Argonne,  IL  60439. 

56.  T.  W.  Armstrong,  Science  Applications,  Inc.,  P.O.  Box  2807,  La  Jolla, 
CA  92038. 

57.  Miguel  Awschalom,  National  Accelerator  Laboratory,  P.  0.  Box  500, 
Batavia,  IL  60510. 


58.  D.  Berley,  National  Science  Foundation,  Washington,  D.C.  20550. 


59.  Larry  W.  Brackenbush,  Battelle  PNL,  P.O.  Box  999,  Richland,  WA 
99352. 

60.  Dr.  J.  E.  Brau,  Dept,  of  Physics,  University  of  Tennessee,  Knoxville 
TN  37916. 

61.  Brookhaven  National  Laboratory,  Attention:  Research  Library,  Upton, 
NY  11973. 

62.  Dr.  Bruce  Brown,  Fermi  National  Accelerator  Laboratory,  P.  0.  Box 
500,  Batavia,  1L  60510. 

63.  Dr.  David  0.  Caldwell,  Department  of  Physics,  University  of 
California  at  Santa  Barbara,  Santa  Barbara,  CA  93106. 

64.  Dr.  Ernest  Coleman,  Office  of  High  Energy  Physics,  U.S.  Department  of 
Energy,  MS  ER-221,  Washington,  D.C.  20545. 

65.  Stanley  B.  Curtis,  Lawrence  Radiation  Laboratory,  Bldg.  29,  Room  213, 
Berkeley,  CA  94720. 

66.  Anna  Di  Cicacco,  Exp.  VAl,  CERN,  CH1211,  Geneva,  Switzerland. 

67.  Herbert  Destaebler,  Stanford  Linear  Accelerator  Center,  Stanford 
University,  Stanford,  CA  94305. 

68.  Robert  T.  Devine,  Defense  Nuclear  Agency,  CDR,  USN,  Washington  D.C. 
20305-1000. 

69.  Dr.  G.  Dietze,  Physikalisch-Technische ,  Bundesallee  100,  D-3300 
Braunschweig,  WESTERN  GERMANY. 

70.  R.  D.  Edge,  Physics  Department,  University  of  South  Carolina, 
Columbia,  SC  29208. 

71.  C.  M.  Eisenhauer,  National  Bureau  of  Standards,  Bldg.  245, 
Gaithersburg,  MD  20899. 

72.  Dr.  R.  Eisenstein,  Department  of  Physics,  University  of  Illinois, 
Urbana,  IL  61801. 

73.  Doug  E.  Engleson,  Defense  Nuclear  Agency,  Armed  Forces  Radiobiology 
Research  Institute,  Bethesda  MD  20814. 

74.  Lt.  Richard  Engstrom,  Air  Force  Weapons  Laboratory,  Kirkland  Air 
Force  Base,  Albuquerque,  NM  87117. 

75.  Dr.  Chris  Fabjan,  CERN,  Geneva  23,  Switzerland. 

76.  Dr.  G.  Feldman,  Stanford  Linear  Accelerator  Center,  Stanford 
University,  Stanford,  CA  94305. 

77.  K.  Goebel,  Health  Physics  Group,  CERN,  1211  Geneva  23,  Switzerland. 

106 


78.  Dr.  D.  Hitlin,  Stanford  Linear  Accelerator  Center,  Stanford 
University,  Stanford,  CA  94305. 

79.  Craig  Jensen,  Dept,  of  Nuclear  Engineering,  University  of  Oklahoma, 
Norman  OK  73019. 

80.  Mr.  Terrence  Jensen,  Dept,  of  Physics  and  Astronomy,  The  Universi  ty 
of  Rochester,  Rochester,  NY  14627. 

81.  Prof.  D.  Lai,  Tata  Institute  of  Fundamental  Research,  National  Centre 
of  the  Government  of  India  for  Nuclear  Science  &  Mathematics,  Homi 
Bhabha  Rd. ,  Bombay  5,  India. 

82.  Lawrence  Livermore  Laboratory,  Technical  Information  Department,  P.0. 
Box  808,  Livermore,  CA  94550. 

83.  Library  for  Nuclear  Science,  Massachusetts  Institute  of  Technology  at 
Middleton,  Middleton,  MA  01949. 

84.  Dr.  V.  S.  Narasimham,  Tata  Institute  of  Fundamental  Research,  Bombay 
400  005,  India. 

85.  W.  R.  Nelson,  Stanford  Linear  Accelerator  Center,  Stanford 
University,  P.  0.  Box  4349,  Stanford,  CA  94305. 

86.  Reran  O'Brien,  Health  and  Safety  Laboratory,  U.S.  Department  of 
Energy,  376  Hudson  Street,  New  York,  NY  10014. 

87.  Dr.  T.  R.  Palfrey,  Jr.,  Department  of  Physics,  Purdue  University, 

West  Lafayette,  IN  47907. 

88.  Dr.  Robert  Palmer,  Brookhaven  National  Laboratory,  Upton,  NY  11973. 

89.  Dr.  C.  W.  Peck,  Department  of  Physics,  California  Institute  of 
Technology,  Pasadena,  CA  91109. 

90.  Dr.  R.  E.  Pevey,  Savannah  River  National  Laboratory,  Aiken  SC  29808. 

91.  Dr.  David  W.  0.  Rogers,  Head,  Ionizing  Radiation  Standards,  Division 
of  Physics,  National  Research  Council  of  Canada,  Ottawa,  Canada  K1A0R6 . 

92.  Dr.  C.  Rubbia,  Lyman  Laboratory,  Harvard  University,  Cambridge,  MA 
02138. 

93.  Dr.  W.  Schmidt,  Institute  of  Experimental  Nuclear  Physics,  University 
of  Karlsruhe,  75  Karlsruhe,  W.  Germany. 

94.  The  Secretary,  Radiation  Group,  Lab  II,  CERN,  1211  Geneva  23, 
Switzerland . 

95.  Dr.  Walter  Selove,  University  of  Pennsylvania,  Department  of  Physics, 
Philadelphia,  PA  19104. 


96.  B.  S.  P.  Shen,  Department  of  Astronomy,  University  of  Pennsylvania, 
Philadelphia,  PA  19104. 

97.  Stanford  Linear  Accelerator  Center,  Attention:  Library,  P.  0.  Box 
4349,  Stanford,  CA  94305. 

98.  Dr.  Alan  Stevens,  Accelerator  Dept.  Bldg.  911,  Brookhaven  National 
Laboratory,  Upton,  NY  11973. 

99.  G.  R.  Stevenson,  Radiation  Protection  Group,  Lab  II,  CERN,  1211 
Geneva  23,  Switzerland. 

100.  Dr.  L.  Sulak,  Department  of  Physics,  University  of  Michigan,  Ann 
Arbor,  MI  48109. 

101.  R.  Tesch,  DESY,  Hamburg,  Notkesteig  1,  W.  Germany. 

102.  Ralph  H.  Thomas,  Lawrence  Berkeley  Laboratory,  Occupational  Health 
Office,  50A-5106,  Engineering  Division,  Berkeley,  CA  94720. 

103.  Dr.  W.  J.  Willis,  CERN,  Geneva  23,  Switzerland. 

104.  Dr.  S.  Yellin,  Stanford  University,  Stanford  Linear  Accelerator 
Center,  P.  0.  Box  4349,  Stanford,  CA  94305. 

105.  G.  B.  Yodh,  University  of  Maryland,  College  Park,  MD  20742. 

106.  Robert  Young,  Director,  Defense  Nuclear  Agency,  ATTN:  STRP,  6801 
Telegraph  Road,  Alexandria,  VA  22310. 

107.  Dr.  B.  Zeitnitz,  Nuclear  Research  Center,  Karlsruhe,  W.  Germany 
93120. 

108.  Gary  H.  Zeman,  Defense  Nuclear  Agency,  CDR,  MSC,  USN,  Armed  Forces 
Radiobiology  Research  Institute,  Bethesda  MD  20814. 

109-343.  Defense  Nuclear  Agency  Transport  Distribution  (AU) . 

344-373.  Technical  Information  Center,  Oak  Ridge  Turnpike,  Oak  Ridge  TN 
37830. 


